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ABSTRACT 


A three—dimensional analysis of the intensity distribution of backward optical 
transition radiation has been performed. The effects of variations in electron energy 
and beam divergence and on material properties such as dielectric permittivities and 
the resultant coherence length upon the angular distribution and polarization of 
optical transition radiation has been investigated. A surprising obeservation 
important to the use of optical transition radiation as a diagnostic tool for high 
energy electron beams is the behavior of the perpendicular component of the 
intensity. In contrast to low energies where the parallel component dominates, at 
electron energies above 200 MeV, the perpendicular component dominates. This 
requires the use of a polarization filter to diagnose particle beam properties at high 
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I. INTRODUCTION 


Transition radiation is produced whenever a uniformly moving charged 
particle passes from one medium into another. In introducing the concept in 1944, 
Frank and Ginsburg [Ref. 1] noted that the intensity, polarization and angular 
distribution of transition radiation are dependent upon the dielectric constants of 
the two media and that transition radiation is not connected with changes in 
velocity of the charged particle. This behavior of transition radiation, particularly 
the angular distribution and polarization dependence on energy makes it a 
promising technology for use as a diagnostic for charged particle beams at low 
energy. 

Since its introduction, many investigations of the nature of transition 
radiation in such diverse environments as the boundary between two media, 
electromagnetic fields. plasmas, inhomogeneous media, and media in which the 
properties vary with time have been published. providing an _ exhaustive 
bibliography. Ginzburg and Tsytovich [Ref. 2], Ter—Mikaelian [Ref. 3], and 
Garibian [Ref. 4] provide brief overviews of some of the established concepts 
resulting from those investigations. Unfortunately, the calculations on the subject 
were always awkward in nature, since the formulas developed were very 
complicated. Many formulas developed independently of each other differed 
considerably from one publication to the next. The situation was exacerbated by 
the lack of specifics or standardization with respect to definitions of key concepts 


and adoption of universal conventions. Wartski [Ref. 5] attempted to rectify the 


problem by returning to fundamental theories, building upon basic hypotheses, and 
emphasizing the conventions adopted by early pioneers in this field. Rule and 
Fiorito et al. [Ref. 6-8] refined and extended Wartski's work to develop analytical 
solutions concerning the components of radiation intensity. Much of the analysis in 
this thesis is based upon those analytical solutions. 

It is necessary to first discuss the physical processes involved in which 
transition radiation occurs. The simplest case to consider is that of a charged 
particle moving at a uniform velocity.which encounters a boundary between two 
media. The media are characterized by their respective dielectric permittivities, q 
and 0: The dielectric permittivities are complex functions of the refractive index 
n and the absorption coefficient « of the medium. The motion of the charged 
particle creates, by definition, an electrical current which in turn generates 
associated electric and magnetic fields. Since the particle is assumed to travel in a 
uniform motion in which it does not experience acceleration, it does not emit 
radiation until encountering the boundary between the media. 

At least three types of radiation are emitted as the particle transits through 
the interface between the media. These are Cherenkov, bremsstrahlung, and 
transition radiation. Transition radiation is closely related to Cherenkov radiation 
but exhibits different properties. Bremsstrahlung radiation is an entirely separate 
mechanism involving particle acceleration. The form of transition radiation may be 
found by examining the fields involved. The electric and magnetic field in each 
medium are described by the solutions to the Maxwell equations for a point charge 
moving at constant velocity. However, the electric and magnetic field components 


in each media fail to satisfy the continuity requirements across the interface. 


The solutions of homogeneous Maxwell equations must be added to those of the 
electric and magnetic fields in each medium to meet the boundary conditions. The 
required radiation fields described by these solutions meeting the boundary 
conditions are the transition radiation. [Ref. 4,9] 

Obviously, finding the solution for the transition radiation fields is quite 
complicated. Ter—Mikaelian [Ref. 3] and Wartski [Ref. 5] provide an excellent 
account of the full derivation. The results of those derivations will be presented and 
analyzed in this thesis so as to promote a better understanding of transition 
radiation and its potential applications for beam diagnostics. 

As stated earlier, the simplest case is that of a charged particle moving at 
uniform velocity encountering a boundary between media. For purposes of this 
Study. a vacuum to medium transition was assumed in which GE 1 and é> Can be 
written simply as ¢. The medium encountered is a thin metallic foil. Transition 
radiation may be observed from the front face of the foil and is called backward 
transition radiation since it is reflected back from the foil. Transition radiation is 
also emitted in the forward direction from the back side of the foil and is thus 
termed forward transition radiation. Observation of forward transition radiation is 
difficult because of its proximity to the beam of charged particles. The backward 
transition radiation. however. may be diverted from the beam by reflecting it at an 
angle. For this reason. the beam of charged particles may be made to encounter the 
metallic foil at an oblique angle of 45 degrees. This has historically been the 
preferred angle of incidence because, since the angle of reflection equals the angle of 


incidence, it allows the observation equipment to be set up 90 degrees from the 


beam axis along the axis of reflection. This configuration permits the study of 
backward transition radiation apart from the beam while keeping calculations 
relatively simple. 

The mean angle at which optical transition radiation appears is the angle of 
specular reflection. If the beam has an angular divergence then the charged 
particles have different angles of incidence upon the transition radiation foil. If the 
angle of a particular particle is a@ with respect to the beam axis, then the optical 
transition radiation will appear at an angle a to the axis of specular reflection. 
Figure | demonstrates the definition of the particle angle a with respect to the 


beam axis, angle of specular reflection W, and angle of observation 0. 


SSS ee >= BEAM AXIS 


tT INTERFACE 
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Figure 1. Definition of the particle angle a with respect to the beam axis, angle of 
specular reflection W, and the observation angle @ with respect to the angle of spec 
reflection for a particle encountering a boundary at an oblique angle. [Ref. 5,9]. 


The interest here is the intensity of the radiation at any point in the plane of 
observation as a function of the observation angle §. The angular information of 
transition radiation obtained from this is valuable as a potential diagnostic for 
charged particle beams. The TR3 program used in this study was based upon the 
analytical solution of this single foil case developed by Rule and Fiorito [Ref. 6-8] 
for this purpose. This program extends Rule and Fiorito's solution into three 
dimensions, thus using the advantage of three-dimensional graphics to aid in 
interpreting the results. 

Recall that the motion of the charged particle induces electromagnetic fields in 
the medium and that radiation is emitted when the particle encounters a boundary 
between two media of dissimilar dielectric constants. When that transition from 
one medium to another occurs, the electromagnetic fields induced by the particle 
motion as well as the fields induced by the transition must adjust to the change in 
dielectric permittivity. The distance over which this adjustment occurs is known as 
the coherence length. (It is also sometimes referred to as the formation length). 
The definition of this length is a bit arbitrary, but is usually chosen so that the 
radiation emitted at the interface adds coherently with the fields generated by the 
particle. It has been found that a phase difference of one radian is the condition for 
the rays to become coherent. The fields at much greater distances may be 
considered to be pure radiation fields. 

Coherence length is related to both the particle velocity and to the phase 
difference between the particle's fields and the radiation fields from the transition. 


The coherence length in the medium can be expressed [Ref. 5] as 
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where (£ is the ratio of the particle velocity to the speed of light, w is the angular 
frequency ( =27c/A), @ is the angle of observation, and <A is the wavelength of the 
radiation in the medium. The coherence length in a vacuum may be found by 
setting /e=1. 

Coherence length takes on a special significance when considering a particle 
encountering more than one boundary as with the two—foil interferometer developed 
by Wartski [Ref. 5]. This interferometer consists of two parallel foils set at an angle 
of 45 degrees to the trajectory of the charged particle beam. Particles traversing the 
foils emit radiation from the back of the first foil that is reflected off the front 
surface of the second foil. Transition radiation is also produced by the transit 
through the second foil. The radiation emitted from the front of the second foil 
forms an interference pattern with the radiation from the back of the first foil 
(Figure 2). The interference pattern is centered around the axis of specular 
reflection. The phase difference between the transition radiation emitted from the 


first and second foils is 


) = (2aL/AB)(1-fcos8) = L/L, (1.2) 


where L is the separation between the foils and L., is the formation length in a 
vacuum. This phase difference should be an integer multiple of one radian for the 
transition radiation from the two foils to add coherently. That is, the separation 


between the foils should be an integer multiple of the coherence length in a vacuum. 
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Figure 2. Diagram of the Wartski's two—foil interferometer. [Ref. 5,9] 


Rule and Fiorito [Ref 6-8] developed an analytical solution to the transition 
radiation emitted froin Wartski's two—foil interferometer [Ref. 5]. The COHER3 
program used in this study extended that solution into three dimensions for use with 
three—dimensional graphics to aid the analysis. 

Most current theories on transition radiation were built on the simplification 
of assuming a point value for both energy and frequency in order to evaluate 
complex integrals. In reality, however, energy and frequency extend over a finite 
range of values. The assumption may be responsible for some slight differences 


between theoretical and actual data. The capability of examining the effects of a 


finite bandwidh of wavelength (and thus frequency) was built into the COHER3 
program. This thesis provides one of the few investigations on finite frequency 
bandwidth effects. 

The TR3 and COHER3 programs created an output compatible with the 
Three—dimensional graphics program, SON OF SURF3D. This graphics program 
was written by Don Gilbert of Dogstar Software as an extension of SON OF SURF. 
At the time of this writing, SON OF SURFS83D was still under development and not 
yet ready for public release. This thesis was one of the first thorough testings of the 
SON OF SURF3D graphics program. As such, it may play a key role in 
determining whether the SON OF SURF3D program will be fully developed for 
commercial use. Use or discussion in this thesis of any comercially availabe product 
does not constitute endorsement. 

The purpose of this thesis was to promote a better understanding of transition 
radiation based upon the theoretical developments presented in references 1—9. It 
provides an investigation of the effects of finite bandwidth on radiation intensity 
distribution. The programs and graphics package used in this thesis could become 
powerful research tools in the development of TR as a diagnostic for charged 
particle beams. The behavior of the polarization and angular distribution of 
transition intensity in response to changes in specific parameters demonstrated in 
this thesis should provide a guide for further development of transition radiation 
theory for this purpose. By building upon basic concepts and examining the nature 
of the complex formulas involved, with the aid of three—dimensional graphics, this 
thesis proposes to clarify the concepts of transition radiation so that even the 


laymen could obtain a basic understanding. 


Chapter II descibes the analytical solutions of Rule and Fiorito for transition 
radiation from the single foil and two—foil interferometer cases {Ref. 6—9] and the 
TR3 and COHER3 programs used in this study which extended those solutions to 
three dimensions. The chapter closes with a brief description of the three 
dimensional graphics program, SON OF SURF3D, used in conjunction with these 
programs. 

Chapter III contains the analysis of the effects of energy, beam divergence, and 
dielectric properties on the distribution of transition radiation intensity for the 
single foil case. An investigation into the effects of coherence length and optical 
frequency bandwidth as well as energy, beam divergence, and dielectric properties 
for the two—foil interferometer case forms the content of Chapter IV. 

The conclusions resulting from the analysis as well as recomendations for 
further study, hardware, and software are discusssed in Chapter V. For 
convenience, a users manual for the single foil and two—foil interferometer programs 
and for the three-dimensional graphics program was developed and included as 
Appendix A. Appendix B contains the listing for the single foil program TR3 and 


the two—foil interferometer program COHER3. 


Il. CALCULATIONAL APPROACH 


The core of the analysis of this study are analytical solutions developed by D. 
W. Rule et al. of the Naval Surface Warfare Center [Ref 6—9] for single foil and for 
two—foil interferometer transition radiation. A series of programs developed by D. 
W. Rule et al. around those analytical solutions were used in the study of transition 
radiation characteristics to produce theoretical two—dimensional profiles of relative 
intensity [Ref. 7—9]. 

Rule's programs were modified to extend the analysis to three dimensions and 
to add flexibility in setting data parameters during run time. The resulting 
program for the single foil case is called TR3 and an interferometer program which 
took into account coherence length and frequency bandwidth is called COHERS. 
Both programs were tailored specifically to be compatible with the 
three—dimensional graphics program, SON OF SURF3D, written by Don Gilbert of 
Dogstar Software [Ref. 10}. SON OF SURF3D was in the developmental stage, not 
yet ready for public release as an improved version of SON OF SURF, extending the 
capabilities of that program to plot three—dimensional surfaces from data files. A 
description of SON OF SURFS3D follows the discussion of the theoretical basis of the 
TR3 and COHERS programs. 


A. SINGLE FOIL 
The program TR3 was a modification of TRADS1, which was one of a series of 
single—foil transition radiation programs written by D. W. Rule [Ref. 8]. The 


TRADS1 program calculated a two—dimensional intensity distribution as a function 
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of angle measured from the angle of specular reflection. The angle of specular 
reflection was chosen to be forty—five degrees as measured from the target plane and 
ninety degrees from the beam axis (Figure 1) for the reasons discussed in chapter 
one. The program assumes a Gaussian distribution of beam divergence angle and a 
vacuum to medium transition. The medium is characterized by a dielectric 
constant which is in general complex. All calculations were based on energy and the 
incremental values of the observation angle (6), measured from the angle of 
specular reflection. The TR3 program extended the analysis of the TRADS1 series 
of programs to three—dimensions in conjunction with the three—dimensional 
graphics program SON OF SURFSD. 

The values of the parallel and perpendicular components of intensity as well as 
the total intensity were based on an analytical approach developed by D. W. Rule 
and associates [Ref. 7,8] in accordance with the theoretical treatment of transition 
radiation by Wartski [Ref. 5]. Using the small angle approximation of cos @ and sin 
). Rule [Ref. 7] showed that, for small angles, the intensity per unit frequency and 


solid angle could be expressed 
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where ( is the velocity of the particle in units of the speed of light, B is the 


component of the velocity that is perpendicular to the observation plane, 7 is the 


1] 


dimensionless Lorentz factor for energy, and cI and rare respectively the parallel 
and perpendicular Fresnel reflection coefficients. The Lorentz factor y and the 


velocity £ are defined 








and 


The parallel and perpendicular reflection coefficients rl and rare calculated 
in terms of the the observation angle @ and dielectric constants of the media. As 
discussed earlier, a vacuum to medium transition was asssumed in which the 
dielectric constant € of the vacuum is equal to one, and the dielectric constant of 


the medium may be written as ¢. Then ri and rare given by 


_ ecos? — ee 
"| > ; 
Ecos = ee 


ma 
a 
cos? + | e-sin2g 


Equations (2.1) and (2.2) contain only the two highest terms in the Lorentz 
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factor y [Ref. 8]. Note that the odd—powered factor of @ in I) gives rise to an 


Ve 


asymetric pattern as the angle changes sign. The TR3 program preserves this 
asymetric behavior by reversing the sign of the angles represented in the first 32 
columns in the 64 by 64 array used to calculate and organize the data. 

Rule's analytical calculations took into account the effect of beam divergence 
on the transition radiation intensity pattern. This was done by folding Gaussian 


distributions for the projected angles into equations (2.1) and (2.2), resulting in 
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where Z is defined as (77 +i8)/ 20, 7, as the rms beam angle of divergence 
projected into the observation plane, and x as the rms angle of beam divergence 
perpendicular to this plane and containing the beam axis. Rule ignored the second 
term of equation (2.2) for equation (2.8) and it was also assumed that the rms beam 
divergence was small, i.e. ee The function W(Z) is given in terms of the 


error function (Z) in the non-standard form 


W(Z) = [1-4(Z)Jexp(Z”). (2.9) 


The whole of the single foil program, TR3, was built around the analytic 


solutions embodied by equations (2.7—2.9). Total intensity was defined simply as 
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the sum of the perpendicular and parallel components. All intensities were 
calculated in terms of charge squared divided by the speed of light per unit 
frequency and solid angle and were normalized for plotting purposes. Actual values 
of intensities as defined here were stored in a file labelled RADOUT.DAT. Note 
that if the beam divergence is symmetric such that TSO then the total 
divergence angle is ¥20. This case was built into the single foil program. 

Rule's analytic solutions were extended into three dimensions by calculating 
the parallel and perpendicular components of radiation intensity and the total 
radiation intensity over an angular distribution centered around the angle of 
specular reflection. This angular range was broken down into a 64 by 64 array in 
which the intensities were calculated for each point in the array. The parallel 
component, perpendicular component, or total intensity for each point may be 
stored in the graphics output file for plotting. Each point in the array represents an 
angle and position measured from the axis of specular reflection. The axis was set 
at the center of the array in the thirty—second column of the thirty—second row of 
the array. In terms of data storage, the rows of the array correspond to records and 
the columns correspond to field elements. Each of the 64 records contain 64 field 
elements. Since the program takes several minutes to calculate the data, progress is 


indicated by a record number. 


B. INTERFEROMETER 

The interferometer program, COHER3, was based on the TEM series of 
transition radiation programs developed by D.W. Rule et al. that culminated in the 
two—dimensional program COHERI. The COHERI program was dveloped around 


Rule's analytical solutions for transition radiation in the case of Wartski's two—foil 
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interferometer. A finite optical frequencycy bandwidth was also built in, allowing 
for a study of the affects of finite bandwidth on transition radiation intensity 
distribution. The COHER3 program simply extended COHER1 to three 
dimensions, modifying it so as to allow flexibility in altering data parameters during 
run time, and tailoring the output to be compatible with the three—dimensional 
graphics program, SON OF SURF3D 

At the time of this writing, Rule and Fiorito had not published a full 
discussion of their analytical solution for the two—foil interferometer case. However, 
a discussion of the development of the analytical code for parallel component of 
intensity was published [Ref. 7] and formed the basis of the presentation here. The 
perpendicular component follows a similar argument in development. 

Rule and Fiorito's analysis began with Wartski's [Ref. 5] development of the 
two-foil interferometer. Wartski showed that the parallel component of TR 
intensity per unit frequency and solid angle in this case could be written 
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where w represents frequency, 2 represents a solid angle, w represents the angle oi 
specular reflection, § represents the observation angle as measured from 7, and 
F(w,6,w) represents a Fresnel coefficient which for a vacuum to medium transition 1s 


defined 
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The exponential term is the interference term for the transition radiation 
patterns differing in phase by @¢. This phase difference is defined in equation (1.2) 
as the ratio of the interfoil spacing to the coherence length. For relativistic 
particles, with a Gaussian distribution of beam angles folded in and averaged over 
particle angles with respect to the beam a, the intensity per unit frequency and 


solid angle may be written 
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The basis of the COHER3 program is the analytical code developed by Rule 
and Fiorito to solve both equation (2.12) and the integral for the perpendicular 


component of intensity as well. The result for each component may be written 
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C. THE GRAPHICS PROGRAM 

The three—dimensional graphics program used in conjuction with the TR3 and 
COHER3 programs for the analytical study of transition radiation was SON OF 
SURF3D [Ref. 10]. This program was a derivative of a program written by Don 
Gilmore of Dogstar Software to explore the graphics capabilities of TurboHalo 3.0 
[Ref. 11]. SON OF SURFS3D is able to plot three—dimensional surfaces from data 
files containing three—dimensional point information of up to 100 by 100 arrays. 

The data points stored in the output file for the graphics program SON OF 
SURF3D are automatically standardized relative to the axes according to the 
following formula [Ref. 10] 


Z' = (2-2 nin 4 cen) . (2.15) 


max min 


The TR3 and COHER3 programs store normalized real values of the parallel 
component. perpendicular component, or total intensity in an output file named by 
the user. For purposes of the analysis, the z-maximum value was set at 1.0 and 
both the z—minimum and z—center values were set at 0.0. Doing so sets the axis 
markers at the bottom of the data set in the three—dimensional plots and shows the 
intensties relative to a maximum intensity. An option was provided to normalize 


all intensities to the maximum intensity at 100 MeV. 
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III SINGLE FOIL TRANSITION RADIATION 


The TR3 program was run many times while changing one parameter each 
time in order to study the effect of that parameter on the distribution of 
transition radiation intensity. The effects of energy, beam divergence, 


scattering foil, and dielectric constant foil were studied in this manner. 


A. ENERGY 

The single foil case for transition radiation was calculated for energies of 1 
GeV, 500 MeV, and for 10 MeV to 100 MeV in increments of 10 MeV. The total 
intensity distribution as a function of angle was calculated over an angular 
range of —.05 to +.05 radians with respect to the angle of specular reflection. 
Several observations were made from the resulting data (Tables I,IJ) and 
associated three—dimensional plots (Figures 3-7). The data for Figures 3-7 
were normalized to the maximum value of intensity for the energy plotted. 
Therefore the maximum relative value is one (1.0), which is represented by the top 
of the z—axis. The units on the z—axis are thus the dimensionless ratio of intensity 
to maximum intensity. The angular distribution was over 0.05 radians as measured 
from the z—axis in either direction along the x and y axes. The endpoints on the x 
and y axes coincide with .0625 radians. 

The most obvious observation was that the shape of the plots varied 
greatly from one energy to the next when using the fixed parameters described 
above. However, the plots appear nearly identical when plotting the data over a 


relative scale such as a multiple of the predicted angle of maximum intensty. 
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Such is the case for Figures 8-10 in which the data was calculated over an angular 
distribution of five times the predicted angle of maximum intensity for the energy 
plotted. i.e. | @| < 5/7. 

At low energy the plot was a gently curved cone shape centered around 
the z—axis with the opening extending almost the full angular width of the 
plot. The angular separation measured from the angle of specular reflection 
(z-axis) was defined as theta (Figure 1 of chapter I). The angle between the 
z-axis and the sides of the cone form a particular value of theta. As the 
energy was increased, the cone became more sharply defined and the angle 
steadily decreased. hWartski [Ref. 5] showed this phenomenon was predicted in 
theory from the equation for transition radiation intensity per unit frequency 


and solid angle in this case for a vacuum to medium transition. 
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where q is the charge on the particle, ¢ is the dielectric constant of the medium, 


and the fresnel coefficients ry and ty are defined as 
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Figure 3. Single foil transition radiation at 10 Mev. Note the overall cone shape of 
the plot. The sides of the cone describe the angle at which the maximum intensity 
occurs relative to the axis of specular reflection (Z—axis). The X—Y plane depicts 
the angular distribution in radians. The data extends to from —.05 to +.05 radians 
on each axis. The endpoints on the X and Y axes correspond to +.0625 radians. 
The units measured along the Z—axis are the dimensionless ratio of intensity to 


maximum intensity. The top of the Z-axis corresponds to the maximum value of 
the data. 
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Figure 4. Single foil transition radiation at 40 Mev. The cone shape still exists in 
the center of the plot, but the angle describing the position of maximum intensity is 
considerably smaller. The X—-Y plane depicts the angular distribution in radians. 
The data extends to from —.05 to +.05 radians on each axis. The endpoints on the 
X and Y axes correspond to +.0625 radians. The units measured along the Z—axis 
are the dimensionless ratio of intensity to maximum intensity. The top of the 
Z—axis corresponds to the maximum value of the data. 





Figure 5. Single foil transition radiation at 70 Mev. The original shape of a cone is 
nearly lost in the center of the plot. The angle describing the position of maximum 
intensity is again much smaller. The X—Y plane depicts the angular distribution in 
radians. The data extends to from —.05 to +.05 radians on each axis. The 
endpoints on the X and Y axes correspond to +.0625 radians. The units measured 
along the Z-axis are the dimensionless ratio of intensity to maximum intensity. 
The top of the Z—axis corresponds to the maximum value of the data. 
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Figure 6. Single foil transition radiation at 100 MeV. The only remaining evidence 
of the original cone appears to be the rim occurring at the top of the plot around the 
Z—axis. The apparent asymmetry of the rim is due to the plotting algorithm and is 
not a physical reality. The X—Y plane depicts the angular distribution in radians. 
The data extends to from —.05 to +.05 radians on each axis. The endpoints on the 
X and Y axes correspond to +.0625 radians. The units measured along the Z—axis 
are the dimensionless ratio of intensity to maximum intensity. The top of the 
Z—axis corresponds to the maximum value of the data. 
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ligure 7. Single foil transition radiation at 1 GeV. The data appears to converge to 
a point below the top of the Z—axis which should represent the maximum value oof 
the data. In fact, the maximum value occured on the axis of specular reflection 
represented by the Z—axis. Furthermore, maximum value of the data occurred as a 
point which is at the top of the Z-axis. The X~-Y plane depicts the angular 
distribution in radians. The data extends to from —.05 to +.05 radians on each axis. 
The endpoints on the X and Y axes correspond to +.0625 radians. The units 
measured along the Z—axis are the dimensionless ratio of intensity to maximum 
intensity. The top of the Z—axis corresponds to the maximum value of the data. 
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Figure 8. Single foil transition radiation at 30 Mev. Compare with the plot for 
single foil transition radiation at 40 MeV (Figure 4). In this case data extends to 
five times the predicted value of maximum intensity, 1/y¥ in either direction on the 
X and Y—axes. The Z—axis still measures the relative intensity of the data with the 
top of the Z—axis corresponding to the maximum intensity. 
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Figure 9. Single foil transition radiation at 60 Mev. The plot is nearly identical to 
the plot for 30 MeV (Figure 8), which is plotted on the same relative scale. That. is, 
the data extends to five times the predicted value of maximum intensity, Lain 
either direction on the X and Y—axes. The Z—axis measures the relative intensity of 
the data with the top of the Z—axis corresponding to the maximum intensity. 
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Figure 10. Single foil transition radiation at 90 Mev. Compare with the plot for 30 
and 60 MeV (Tigures 8,9), measured on the same relative scale. The plots appear to 
be identical. The shape of the single foil transition radiation pattern is independent 
to the energy. Intensity increases with energy increases while the angle at which the 
maximum intensity occurs decreases. The data extends to five times the predicted 
value of maximum intensity, 1/y in either direction on the X and Y—axes. The 
Z—axis still measures the relative intensity of the data with the top of the Z—axis 
corresponding to the maximum intensity. 
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The fresnel coefficients must satisfy the continuity relationship: 


For metallic media in the visible region, |e|>>1, and equation 


(3.1) can be rewritten 
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A second approximation was made by considering the angle theta 
sufficiently small such that sin’ 6 and cos“ 0 could be replaced by 6? and 


? 
(1-0/2), respectively, leaving 
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Again, the preceeding equations were for the specific case at hand of a 
vacuum to medium transition. The factor ry’ the reflectance of the 
medium, can be ignored for a medium to vacuum transition. Equation (3.7) 
showed the relationship of radiation intensity to energy and the angle 
theta. The nature of this relationship was shown by differentiating 


equation (3.7) with respect to theta: 


| (3.8) 
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This equation demonstrated that the extrema occur for the theta equal 
zero direction and for theta equal 1/7. Using these values for theta in equation 
(3.7) showed that the radiation intensity is zero in the theta equal zero direction 
and is at a maximum for the theta equal 1/ydirection. This angle of maximum 
radiation intensity decreased as energy increased. A look at Figures 3-7 
Showed that the relationship was qualitatively true. A close look at the 
data (Table I). in which the angle measured where maximum intensity actually 
occurred was compared with 1/y demonstrated excellent correlation. 

The energy. the Lorentz factor 7, predicted angle of maximum intensity 
1/7. and measured angle of maximum intensity 8, for a beam divergence of .001 
radians were tabulated in Table JI. The maximum value of the total intensity 
ITOT, maximum value of parallel intensity IPAR, maximum value of 
perpendicular intensity IPERP, the value of parallel intensity at maximum total 
intensity IPAR(O, er and the value of perpendicular intensity at maximum total 


intensity IPERP(@_), were tabulated in Table II for comparison. Intensity 


was calculated per unit frequency and solid angle in terms of the square of 
the charge divided by the speed of light. Energy was listed in units of MeV, 
angles were listed in units of radians. Additionally, parallel, perpendicular, 
and total intensities at the angle of maximum intensity as well as the maximum 
of the perpendicular intensity versus energy were plotted in Figure 11. The 
plot revealed that total and parallel intensity increased as a power 


function of energy. 


TABLE 1.ACTUAL VS. ANGLE OF MAXIMUM INTENSITY 


Energy ¥ ia Actual Yes 
(MeV) (radians) (radians) 
10 201 0486 0487 
20 40.14 0249 0250 
30 59.71 .0168 .0169 
40 79.28 .0126 .0126 
90 98.85 .0101 .0102 
60 118.4 .0084 0087 
70 137.0 .0073 0072 
SO amine’ .0064 .0065 
90 Lid .0057 0058 
100 196.7 0051 0051 
500 979.5 .0010 0000 


It was interesting to note that the maximum intensity at 


.5 GeV and 1.0 


GeV occurred at the theta equal zero direction. A look at the data in Table II 
provides a clue as to why this occurred. As energy increased, the perpendicular 
component of intensity provided a greater proportion of the total intensity. It was 
found that the perpendicular component of intensity began to dominate over the 
parallel component at energies above 200 MeV. At .5 GeV and 1.0 GeV, the 


perpendicular component of intensity dominated over the parallel component of 
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Figure 11. A plot of the components of intensity versus energy. It is apparent that 
the perpendular component at the angle where the maximum intensity occurs 1S 
negligible at lower energies. However, the maximum value of the perpendicular 
component can not be ignored. This maximum value occurs on the axis of specular 
reflection, and thus tends to fill up the cone shape seen in Figures 3—o. 
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intensity in contributing to the total intensity. 
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TABLE ILINTENSITIES 


IPAR 


ctf 


9.567 
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The components of radiation 


intensity behave quite differently from each other and contribute distinct profiles to 


the total intensity distribution. Wartski [Ref. 5], (Figure 12) showed that the 


perpendicular and parallel components of radiation intensity could be 


KEltten: 
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, are as depicted in Figure 12 and By and B. are the 


projections of the %/c vector on the plane of observation. Wartski [Ref. 5] 
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Preure 12. The the unit nea n", and B vectors in the plane of observation as 
depicted in Wartski [Ref. 5]. 


showed that for ultra—relativistic particles, where it is assumed that Je>>1 
in a vacuum to medium transition, the second term in equations (3.9) and 


: + : 
(3.10) dominates as /-n becomes close to one. In this case 
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The maximum value of I does not occur at the same position as the 
maximum of total intensity except at very high energies where the perpendicular 
component of intensity becomes dominate and thus determines the location of the 
maximum total intensity. It was found that the energy where this occurs is in the 
neighborhood of about 200 MeV. A look at the perpendicular reflection coefficient 
reveals that the maximum of the perpendicular component of radiation intensity 
occurs when theta (@),the angular separation from the observation angle (z—axis) is 


ZETO 


| a Oe = 
i 0% ¥ Case il” (en ee . (3.13) 
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For non-relativistic particles, the perpendicular component of the 
radiation intensity could be shown to be essentially zero while the parallel 


component of radiation intensity could be written 
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The £6, in equation (3.14) indicates that the radiation produced by a 


non-relativistic particle for an oblique angle is the same as that for 
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normal incidence at a velocity of B, rather than £@. It was clear from 
equation (3.14) that the parallel component of radiation intensity for a 
non-relativistic particle is proportional to the square of the energy as is 
the total intensity. 

In the ultra—relativistic energy region, fT remains proportional to 
the square of the energy. Equation (3.12) indicates that I is also proportional to 
the square of the energy for lower energies. The contribution of I to the total 
energy is negligible at lower energies but becomes quite significant as energy 
increases. and increases asymptotically as B-n approaches one. Figures 13—15 depict 
the perpendicular component of radiation intensity at 100, 250, and 500 MeV. 

Analysis of equation (3.7) revealed that total intensity was 
proportional to the square of the energy. Figures 16-18 dramatically showed 
the dependence of radiation intensity on _ energy. The intensity 
distribution for these figures was scaled to the maximum intensity value of 
100 MeV. The energies for these figures were 50 MeY to 90 MeV in increments 
of 20 Mel. 

Finally. a closer look at Figures 8-10 reveals that the plots aren't quite 
identical in appearance. In each case the intensity distribution was calculated over 
an angular distribution of five times the predicted angle of maximum intensity of 
1/7 for the energies of 30, 60. and 90 MeV. Closer inspection revealed that the apex 
of the cone decreased in depth as energy increased. Obviously, the apex of the cone 
was not zero intensity at higher energies. This was due to the contribution of the 


perpendicular component of intensity to the total intensity distribution. 
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Figure 16. Single foil transition radiation intensity at 50 MeV. The data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y—axes. The Z—axis measures the intensity of the data relative to the 
maximum intensity obtained for 100 MeV using the default values of the program. 
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Figure 17. Single foil transition radiation intensity at 70 MeV. The data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y—axes. The Z—axis measures the intensity of the data relative to the 
maximum intensity obtained for 100 MeV using the default values of the program. 
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Figure 18. Single foil transition radiation intensity at 90 MeV. The data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y—axes. The Z—axis measures the intensity of the data relative to the 
maximum intensity obtained for 100 MeV using the default values of the program. 
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TABLE III, MULTIPLES OF 1/7 


Energy 1 IL 5/7 ee 
(MeV) (radians) (radians) (radians) 
10 20.57 .0486 .2431 0487 
30 99.71 .0168 .0837 .0168 
60 118.4 .0084 .0422 0087 
90 ive 0057 .0282 .0058 


In summary, the angle of maximum intensity, total intensity, perpendicular 
componenet of intensity, and the parallel component of intensity were found to be 
very sensitive to changes in energy. The perpendicular component of intensity was 
found to be rather insignificant at energies below 50 MeV but quite dominate at 
energies above 200 MeV. The angle of maximum intensity, total intensity, and 
parallel component of radiation intensity would be the most useful parameters for 
diagnosic purposes of beam quality at lower energies. However the angular 
information is lost above 200 MeV when the perpendicular component of intensity 


becomes dominate. 


B. BEAM DIVERGENCE 

The effect of beam divergence upon the transition radiation intensity 
distribution may be studied by assuming a Gaussian distribution of divergence 
angles. The Gaussian distribution of angles may be approximated by taking the 
root mean square average value and applying it to the transition radiation intensity 


equations (2.7) and (2.8) discussed in Chapter II. 
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Recall that when writing the code for the single foil case, Rule [Ref. 7] 

programmed in the assumption that o. = o. = o such that the total divergence 


x y 
angle is Yo. In this case, equations (2.7) and (2.8) may be rearranged to bring out 


the dependence on beam divergence more clearly. 
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It is now clear that the parallel component of intensity should decrease with 
an increase in beam divergence. On the other hand, any increase in beam 
divergence should bring about a corresponding increase for the perpendicular 
component of intensity. Recall that the maximum intensity of the perpendicular 
component occurs on the axis of specular reflection. The maximum of the parallel 


component, however, occurs away from that axis at an angle equal to the inverse of 
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the Lorentz factor. Suppose the beam divergence would be permitted to become 
large enough for the perpendicular component to become comparable with, or even 
dominant over the parallel component. The angle at which the maximum intensity 
occurs would shift, eventually ocurring on the axis of specular reflection whenever 
the perpendicular component dominates over the parallel component. The angular 
information provided by the parallel component would then be lost. Both 
components of intensity increase as a power function of energy so beam divergence 
effects should be more pronounced at higher energies. 

The data in Table IV and in Figures 19—28 bear out the correctness of this 
analysis. ‘Table IV contains a comparison of the total, parallel, and perpendicular 
intensities for various energies and rms beam divergence angles. Energy is listed in 
MeV, beam divergence, labelled Beam Div., is listed in radians. The intensities are 
per unit frequency and solid angle in terms of charge squared divided by the speed 
of light. The values for the parallel and perpendicular components of intensity are 
the maximum values obtained for that component. Since the maxima of the 
components of intensity do not occur in the same position, the sum of the values 
listed will not agree with the value of the total intensity. The angle o. at which 
the maximum intensity occurs 1s in radians. It is included here to study the 
behavior of this parameter as the parallel and perpendicular components of intensity 
adjust roles to changes in beam divergence. Note that an increase in beam 
divergence 1s accompanied by a corresponding increase in the total and perpendicuar 
intensities and a decrease in the parallel component of intensity. Whenever the 


beam divergence was large enough for the perpendicular component to dominate 
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over the parallel component, the maximum intensity occurred on the axis. This can 
be seen in Table IV in that whenever the value of IPERP was larger than that of 
IPAR, the value of ae was Zero. 

The data for Figures 19-28 was normalized to the maximum intensity for the 
energy plotted. The top of the z—axis corresponds to the value of the maximum 
intensty. The horizontal width of the data in the plots of Figures 19-26 was fixed 
at 0.05 radians as measured in either direction along the axis from the z-axis. The 
ends of the X and Y axes correspond to a value of .0625 radians. The horizontal 
width of the data in the plots of Figures 27 and 28 was set at a value of five times 
the predicted angle of maximum intensity. This option was chosen to preserve 
clarity at high energies. 

The effects of beam divergence on the intensity distribution are clearly seen in 
Figures 19—28. Figures 19—23 depict the effects of increasing beam divergence at 10 
MeV. A general flattening of the intensity distribution seems to occur. However, 
recall that the predicted angle of maximum intensity for this energy is .049 radians 
while the data is calculated out to .05 radians. Thus Figures 19—23 provide a close 
look at the region where the perpendicular component of intensity becomes more 
and more dominant. The perpendicular component of intensity appears to become 
completely dominant when the beam divergence was increased to Just over .0250 
radians. 

Figures 24—26 depict the shape of a gaussian distribution at an energy of 40 
MeV as the beam divergence is increased from 0.0010 to 0.0060 radians. In this case 
the predicted angle of maximum intensity at .0126 radians is much less than the 
.0500 radians over which the data is calculated. Therefore, the effects of beam 


divergence on the parallel, as well as the perpendicular, component of intensity can 
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Figure 19. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 1.000E—3 radians. The data extends to .0o radians in the X—-Y 
plane. The top of the Z—axis corresponds to the maximum intensity of the data. 
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Figure 20. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 1.000E—2 radians. The depth of the cone is much shallower than 
in the previous figure. The data extends to .0o radians in the X—Y plane. The top 
of the Z—axis corresponds to the maximum intensity of the data. 
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Figure 21. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 1.250E—2 radians. Note that the depth of the cone continues to 
decrease with an increase in beam divergence. The data extends to .05 radians in 
the X—Y plane. The top of the Z-axis corresponds to the maxirnum intensity of the 


data: 
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Figure 22. Single foil transition radiation pattern at 10 MeV for an rms beain 
divergence angle of 2.000E—2 radians. What's left of the cone appears to be little 
more than a dimple. Note that the angle at which the maximum intensity occurs fo 
an energy of 10 MeV is .049 radians while the data presented here extends to .05 
radians in the X-Y plane. The region depicted then coincides with the region in 
which the perpendicular component of intensity becomes more and more dominant 
over the parallel component. The top of the Z—axis corresponds to the maximum 


intensity of the data. 
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Figure 23. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 2.500E—2 radians. All the angular information of the TR 
pattern has been completely washed out by the effects of beam divergence. The 
data extends to .05 radians in the X—Y plane. The top of the Z-axis corresponds to 
the maximum intensity of the data. 
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Figure 24. Single foil transition radiation pattern at 40 MeV for an rms beam 
divergence angle of 1.000E-3 radians. The maximum intensity occurs at an angle of 
0126 radians while the data presented extends to .05 radians in the X—Y plane. 
The top of the Z-axis corresponds to the maximum intensity of the data. 
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Figure 25. Single foil transition radiation pattern at 40 MeV for an rins beam 
divergence angle of 5.000E—3 radians. Note that the rim around the cone is more 
rounded thatn in the previous figure and that the depth of the cone has significantly 
decreased. The increase in the perpendicular component, which is at a maximum on 
the axis, is responsible for the decreasing depth of the cone. The data extends to .05 
radians in the X—Y plane. The top of the Z—axis corresponds to the maximum 
intensity of the data. 
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Figure 26. Single foil transition radiation pattern at 40 MeV for an rms beam 
divergence angle of 6.000E—3 radians. The angular information has been nearly 
washed out by the beam divergence. The perpendiular component of intensity is 
nearly equal to the parallel component of intensity. The data extends to .05 radians 
in _ X—-Y plane. The top of the Z—axis corresponds to the maximum intensity of 
the data. 





Figure 27. Single foil transition radiation pattern at 100 MeV for an rms beam 
divergence angle of 1.500E-3 radians. The data extends to five times the predicted 
angle of maximum intensity, i.e. 5/7 in the X-Y plane. This option was chosen to 
preserve clarity at this high energy. The top of the Z-axis corresponds to the 
maximum intensity of the data. 
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Figure 28. Single foil transition radiation pattern at 100 MeV for an rms beam 
divergence angle of 3.000E—3 radians. All the angular information has been washed 
out by the effects of beam divergence. Note the apparent gaussian distribution of 
the data. The data extends to +5/y radians in the X—Y plane. The top of the 
Z—axis corresponds to the maximum intensity of the data. 
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TABLE IV, BEAM DIVERGENCE VS INTENSITIES 


Energy Beam Div. OnE IPAR [PE ie 
Di y 2 
(MeV) (radians) ena Sal an] (radians) 
10 1.000E—3 9.571 9.567 1.711E—2 .0487 
10 2.000E—3 Soe | 9.594 6.826E—2 .0485 
10 3.000E—3 CBaie: 9.534 0.1529 0488 
10 4.000E-3 9.574 9.506 0.2703 .0488 
10 5.000E—3 oan 9.471 0.4193 0488 
10 6.000E—3 9.582 9.428 0.5987 0490 
10 7.000E—3 9.588 9.378 0.8069 0490 
10 §.000E—3 9.597 9.322 1.042 .0490 
10 9.000E—3 9.609 9.259 1.304 .0491 
10 1.000E—2 9.625 9.191 1.589 .0490 
10 2.000E—2 OES Se Seas 5.425 0423 
10 2.500E—2 12.00 7.853 1.781 0178 
10 3.000E—2 15.11 1.424 10.29 0.000 
40 2.500E—3 14a Love 22207 .0126 
40 5.000E—3 229 124.9 75.96 .0114 
40 6.000E—3 166.4 119.5 102.4 0082 
40 T.500E—3 214.8 111.9 144.9 0.000 
100 5.000E4 877.8 868.7 Sora 0051 
100 1.500E—3 895.2 §10.2 281.1 0050 
100 2.000E—3 937.1 00 461.7 .0046 
100 3.000E—3 Neo Ju) 690.3 851.6 0.000 
be seen. In Figures 24—26 the parallel component of intensity increasingly 


dominates the distribution at and near the z-axis while the parallel component of 
intensity continues to dominate the intensity distribution farther out from the axis. 
The overall intensity distribution approximates a Gaussian distribution once the 
beam divergence is large enough for the perpendicular component to provide the 
maximum overall intensity on the axis of specular reflection. The point at which 
this occurs for 40 MeV is for an rms beam divergence angle of just over 0.060 

radians. For 100 MeV, this point occured between 0.0020 and 0.0030 radians. 
Figures 27 and 28 display the intensity distribution at 100 MeV for rms beam 


divergence angles of 0.0015 and 0.0030 radians. 
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In summary, increased beam divergence increases the perpendicular 
component of intensity while decreasing the parallel component. The angle of 
maximum intensity decreases as the parallel and perpendicular components adjust 
in response to the increase in beam divergence. This can result in the eventual loss 
of all angular information as the intensity distribution approximates a Gaussian 
distribution. The effect is more pronounced at higher energies and thus a beam 
divergence that may be negligible at low energies may have significant effect at 
higher energies. Beam divergence has a detrimental affect on beam quality, 
especially at higher energies and should be minimized to preserve angular 


information from transition radiation. 


C. DIELECTRIC CONSTANT 

As has been stated before, transition radiation occurs when a charged particle 
in motion encounters a boundary between media. The intensity of the radiation 
emitted is dependent upon the dielectric properties of the media involved. The 
majority of the study of transition radiation here has assumed a vacuum to medium 
transtion. That assumption simplifies the analysis and the formulas involved. Now, 
however, it is more useful to study the equation for the total intensity per unit 
frequency and solid angle in which neither medium is assumed to be a vacuum. 

The classic case, as developed by Frank and Ginsberg [Ref. 1],.1s to assume 
that the particle travelling through a medium which is characterized by adielectric 
constant €, transits an interface to another medium characterized by a dielectric 


constant — €5. The trajectory of the particle is assumed to be normal to the 


Of 


interface. Frank and Ginsberg [Ref. 1] developed expressions for the transition 
radiation intensity in both media. The intensity observed in the first medium per 


unit frequency and solid angle at an angle 0, from the the normal is 


we Ve, sin” 0, cos” 6 2 
a (3.16) 
= oe el 


(1-6 €, +By/eq-e,8in“ 6, ) 
(1-6 ,cos”6, )(1—By €)—€,8in“9,) (€,c080, +Ve, €,—€43in 8, 

The intensity per unit frequency and solid angle observed in the second 
medium at an angle I, to the normal may be found from the above expression by 
substituting —@for @ and interchanging the subscripts 1 and 2. In either case it is 
clear that the intensity observed in the medium is proportional to the square of the 
difference of the dielectric constants of the two media. 

Actually the dielectric constant isn't really a constant at all. It is really a 
slowly—varying function of wavelength that may be considered a constant when 
considering short ranges of wavelength. The dielectric constant is complex, 
consisting of both a real and imaginary part that are related to the refractive index 
n and the absorption coeffiecient « of the medium. Letting ¢1 represent the real 
part of the dielectric constant ¢€ and ¢€2 the imaginary part, the following relations 


define the dielectric constant; 
c= el+ie, (3.17) 


es) 
dl=n" = a0 (3.18) 
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Values for the refractive index and absorption coefficient at discrete 
wavelengths for some metals were found in the Handbook of Chemistry and Physics 
(Ref. 12]. The available information on the dielectric constants of the metals 
limited the data analysis on this subject to three significant figures. TABLE V lists 
the optical wavelength, refractive index n, absorption coefficient «, real part of 
the dielectric constant ¢€1, and the imaginary part of the dielectric constant ¢€2, for 
aluminum Al, gold Au, and silver Ag. Note that the value listed the real part of 
the dielectric constant is the absolute value of the result of equation (3.18). 

Some of the results obtained by this method required verification. The 
formulas listed in equations (3.18) and (3.19) were also found in Wartski [Ref. 5], 
and in Bennett and Bennett [Ref. 13] which also provided a table of optical 
constants at infrared wavelengths for gold, silver, and aluminum. The dielectric 
constants calculated from that table for the infrared wavelengths appeared to be 
consistent with those obtained in Table V. 

Table VI provides the results of applying the dielectric constants listed in 
Table V to the single foil case at an energy of 100 MeV. Recall that the total 
intensity should be proportional to the square of the modulus of the difference of the 
dielectric contants of the two media. For the single foil case, the first media was 
assumed to be a vacuum which is characterized by dielectric constant of one, ie ¢« =1 
+ i0. The total intensity ITOT per unit frequency and solid angle in terms of 
charge squared divided by the speed of light is listed for aluminum Al, gold Au, 
and silver Ag. as a function of wavelength (listed in angstroms), and the square of 


the modulus of the difference between the dielectric constant of the metal and of the 
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TABLE V. DIELECTRIC CONSTANT AT OPTICAL WAVELENGTHS 


Metal Wavelength n K él €2 
(angstroms) 

Al 5590 1.44 3.69 11.5 10.6 

Au 4410 1.18 1.56 1.04 3.70 

Au 5890 0.47 6.02 | 36.0 5.66 

Ag 5000 Osha 71 292 5.88 

Ag 5890 0.18 20.6 424 EZ 


vacuum (6-4) The dielectric constant of the metal for the wavelength listed is 


the same as that listed for the same wavelength in Table V. 


TABLE VI. SINGLE FOIL DIELECTRIC CONSTANT Va) [OTA 
TRANSITION RADIATION INTENSITY AT 100 MEV. 


Metal Wavelength (€,-€y)* oak 
Z 
(angstroms) banal 

Al 5590 138 881 
Au 4410 5.08 654 
Au 5890 1250 $45 
Ag 5000 84700 964 
Ag 5890 179000 970 


The data in Table VI clearly indicates that the total intensity does indeed 
increase whenever the square of the modulus of the difference in dielectric constants 
is increased. Figures 29-31 depict the intensity distributions of aluminum, gold, 
and silver at a wavelength of 5890 angstroms and energy of 100 MeV. The data in 
each plot is normalized to the maximum value, represented by the top of the 
z—axis. The data extends five times the predicted angle of maximum intensity in the 
x-y plane. Since the data in each case is plotted on the same relative scale, Figures 


29—31 appear nearly identical. 
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igure 29. The single foil transition radiation intensity distribution for aluminum 


at 100 MeV at a wavelength of 5890 A. The data extends to +5/y radians in the 
X—-Y plane. The top of the Z-axis corresponds to the maximum intensity of the 
data. 
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Figure 30. ‘The single foil transition radiation intensity distribution for gold at 100 


MeV at a wavelength of 5890 RK. The data extends to +5/y radians in the X—Y 
plane. The top of the Z—axis corresponds to the maximum intensity of the data. 
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Figure 31. The single foil transition radiation intensity distribution for silver at 100 


MeV at a wavelength of 5890 KX. Note that this plot is nearly identical to the plots 
for aluminum and gold (Figures 29, 30). Dielectric permittivity does not affect the 
angular distribution of transition radiation. It does have a significant effect on the 
intensity. The data extends to +5/7 radians in the X-Y plane. The top of the 
Z—axis corresponds to the maximum intensity of the data. 
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IV. COHERENT TRANSITION RADIATION 


The two—foil interferometer program COHER3 was used to study the effects 
of energy, beam divergence, and dielectric constant on coherent transition radiation 
for comparison with the results obtained from the single foil case. The 
interferometer program permits the study of a finite optical bandwidth of 
wavelength and thus also of frequency. In addition, it promotes a_ better 
understanding of the concept of coherence length and its importance in determining 
the separation between foils to obtain coherent addition of transistion radiation 


from successive foils. 


A. COHERENCE LENGTH 

It is necessary to begin the study of the two—foil interferometer with a 
discussion of coherence length. This is due to the fact that the distance between the 
foils of the interferometer in relation to the coherence length, has a profound effect 
upon the resulting interference pattern. Recall from chapter one that the coherence 


length in a medium is defined 


ae ee (4.1) 
27 |1—BYecos8| 


Chapter I also described how the forward transition radiation from the back 
of the first foi! and the backward transition radiation emitted from the front surface 


of the second foil differ in phase by 


64 


b= ($5#)(1-Scos8) = L/L, (4.2) 


where L is the distance between the foils and i is the coherence length in a 
vacuum. Obviously, the interference pattern generated by the interferometer is 
highly dependent upon the phase difference and thus the ratio of the distance 
between the foils and the coherence length. For the radiation patterns to add 
coherently. the distance between the foils should be an integer multiple of the 
coherence length in a vacuum for that energy. Wartski (Ref. 5] showed that the 
coherence length on the axis of specular reflection, i.e. 6 = 0, 1s proportional to the 


square of the energv 


i 
IO 


= 
| 
tw 
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a 
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Table VII displavs the coherence length on the axis of specular reflection at 
various energies for a wavelength of A = 5890 angstroms. This wavelength was 
chosen because it was the only wavelength in the visible spectrum for which the 
optical constants for determining the dielectric constant were available for gold. 
silver, and aluminum. The default value of the two—foil interferometer program for 
the distance between the foils is 1.20 centimeters. 

Figures 32 and 33 demonstrate the impact of coherence length on the 
interference pattern from a two-foil interferometer. For both cases the energy is 30 
MeV and the data in the X—Y plane is spread over +3/y¥ radians from the Z—axis. 
The top of the Z—axis corresponds to the maximum value of the data. In Figure 32 
the distance between the foils is 1.20 centimeters while in Figure 33 the distance is 


equal to the coherence length of .034 centimeters. All angular information available 
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Figure 32. ‘Two-foil interferometer interference pattern at 30 Mev for a clear front 
foil and gold back foil. Distance between foils is 1.20 cm. while coherence length is 


034 cm. at 30 MeV for a wavelength of 5890 A. 
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Figure 33. Two—foil interferometer interference pattern at 30 Mev for a clear front 
foil and gold back foil. Distance between foils is equal to the coherence length of 


034 cm. at 30 MeV for a wavelength of 5890 A. 
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in Figure 33, where the distance between the foils is equal to the coherence length, is 
lost in Figure 32. Not suprisingly the intensities for each case were completely 
different. A numerical comparison of the intensities would not be meaningful since 


one case is coherent and the other is not. 


TABLE VII. COHERENCE LENGTH VS. ENERGY FOR \ = 5890 A 


Energy Coherence Length 
(MeV) (cm) 
10 .004 
20 015 
30 .034 
40 .059 
50 .092 
60 131 
70 aS 
80 VES 
90 .294 
100 363 
150 813 
200 1.44 


B. ENERGY 

The discussion on coherence length demonstrated the effect of energy on the 
considerations for determining foil separation as a function of energy and coherence 
length. Table VIII demonstrates the effect of energy on the angle at which 
maximum intensity occurs, the parallel and perpendicular components of intensity 
and the total intensity for comparison with the effects of the single foil case. The 
intensities listed are for the maximum of each component and therefore should not 
be summed to obtain the maximum total intensity, ITOT. The angle at which the 


maximum intensity occured is much greater than the 1/7 angle predicted for the 
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TABLE VHI. ENERGY EFFECTS ON INTERFEROMETER 
TRANSITION RADIATION 


Energy ey Actual On4 ITOT IPAR IPERP 
Z 2 Z 
(MeV) (radians) (radians) at to | fend 
20 .0249 .0730 31.92 31.92 .O008 
60 0084 .0536 900.3 900.3 1615 
100 .0051 Q277 2046 2046 8868 


single foil case. A comparison with the intensities in Table II of Chapter III 
indicates that the intensities for the interferometer case are much more sensitive to 
increases in energy than for the single foil case. The fact that the intensity of TR 
can be increased through the use of a stack of plates is well known and has proven 
to be a useful tool in studying transition radiation (Ref. 2-4]. Ginsburg and 
Tsytovich [Ref. 1] and Garibian [Ref. 4] showed that a condition of resonance occurs 
when the distance between successive foils is such that the transition radiation fields 
from the foils add coherently. This condition is met when the distance between the 
foils is an integer multiple of the coherence length. If the distance between the foils 
does not meet this criteria, the intensity is then a function of the number of 
boundaries per unit length. For this study the separation distance between the foils 
was set equal to the vacuum coherence length for the energy listed. 

Figures 34—36 depict coherent TR patterns at 60 and 100 MeV. The data is 
spread over an angular distribution of five times the predicted angle of maximum 
intensity for the single foil case. The top of the Z-axis corresponds to the 


maximimum value of the total intensity. 
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Figure 34. Coherent TR interference pattern at 60 MeV for a wavelength of 5890 A. 
separation distance between foils was 1.31 mm. The Z—axis measures relative 
intensity. The top of the Z—axis corresponds to the maximum intensity in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 
Angular information is depicted in the X-Y plane. The data extends to five times 
the predicted angle of maximum intensity (in radians) for the single foil case. 
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Figure 35. Another view of coherent TR interference pattern at 60 MeV for a 


wavelength of 5890 A. Separation distance between foils was 1.31 mm. The Z—axis 
measures relative intensity. The top of the Z—axis corresponds to the maximum 
intensity in units of charge squared divided by the speed of light per unit frequency 
and solid angle. Angular information is depicted in the X—Y plane. The data 
extends to five times the predicted angle of maximum intensity (in radians) for the 
single foil case. 





igure 36. Coherent TR interference pattern at 100 MeV for a wavelength of 5890 


A. Separation distance between foils was 3.63 mm. The Z—axis measures relative 
intensity. The top of the Z-axis corresponds to the maximum intensity in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 
Angular information is depicted in the X—-Y plane. The data extends to five times 
the predicted angle of maximum intensity (in radians) for the single foil case. 
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The characteristics of the surface plot are quite different than that of the 
single foil case. A series of fringes surround the central cone of maximum intensity. 
The positions of these relative maximum and minimum intensities are described by 


Wartski [Ref. 5] in terms of an interference order p which is defined 


p = x4(1-Beos6) = 547-, (4.4) 


V 


where L,, is the coherence length in a vacuum and L is the distance between the 


foils. For the relativistic case, this may be approximated as 


a aD, 


Local minima, gy: occur whenever p is an integer, k. Local maximum 
intensities, te occur whenever p = k+1/2. The order of interference at the center, 
that is on the axis of specular reflection is defined to be 


ae (4.6) 


The angle for which the local maxima and minima intensities occur can be 


described by 


2 n 
Am = “L(PPo) > 4-7) 


for p = k+1/2 and p = k, respectively. 
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The intensites for the coherent interferometer case increase with an increase in 
energy. The intensities for the incoherent case depend upon the number of 
boundaries per unit length as well as upon intensity. The maximum intensity 
occurs at the inner most fringe which is located at a much greater angle than that 
predicted for the single foil case. The locations of maxima and minima are 
described in terms of an order of interference p, wavelength 4A, and distance 
between foils L. The order of interference was shown to be inversely proportional 


to the square of the Lorentz factor 7. 


C. BEAM DIVERGENCE 

The effects of beam divergence were studied for energies of 10, 40, and 100 
MeV for comparison with beam divergence effects at those energies for the single foil 
case. Although an increase in beam divergence decreases the visibility of the 
fringes, the beam divergence required to wipe out all angular information was found 
to be much greater than that required for the single foil case. Coherent transition 
radiation is less susceptible to the effects of beam divergence because of the 
amplified intensities due to resonance. 

Table IX contains a comparison of the total, parallel, and perpendicular 
intensities for a few rms beam divergence angles. In all cases, the rms beam 
divergence angle in the X—Z plane and in the Y—Z plane were set equal to each 
other, i.e. ae In this case the total beam divergence angle is /20, which is the 
value listed in the table in radians. Energy is listed in MeV, the intensities are per 


unit frequency and solid angle in terms of charge squared divided by the speed of 
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light. As in the single foil case, the values listed for each component of intensity are 
the maximum attained for that component and do not necessarily add up to the 


maximum overall intensity, labelled ITOT. 


TABLE IX. BEAM DIVERGENCE VS COHERENT OPTICAL TRANSITI 
RADIATION INTENSITIES AT \ = 5890 & 


Energy Beam Div. Orr IPAR ei RP Actual Ue 
2 2 2 
(MeV) (radians) saa] Coen s 5| (radians) 
10 1.41E—2 Ded Zl 1195 .2604 
10 2.83E—2 alley 20.79 2.685 .2654 
10 4,24E—2 Osi 18.14 5.690 Ana 
10 5.66E—2 ANZ 15.76 9.356 2081 
10 8.49 E—2 26.75 12.50 Louis .O000 
40 1.01E-2 308.2 293.6 2 0688 
40 §.49E—2 300.1 24Se1 51.94 .0695 
40 1.41E-1 2587 69.93 Dott .0000 
100 2.83E—3 es fF 1788 224.9 0277 
100 4,.24E—3 1768 1560 Ae 0274 
100 7.07E—3 1889 1186 1149 .0270 


An increase in beam divergence resulted in a corresponding increase in the 
perpendicular component of intensity and decrease in the parallel component of 
intensity. Also the total intensity tended to decrease until the perpendicular 
component of intensity became comparable to or greater than the parallel 
component. The maximum intensity occurred on the axis of specular reflection 
when the perpendicular component of intensity dominated the _ intensity 
distribution. The physics behind this was explained for the single foil beam 


divergence in Chapter III. 


(is) 


The effect of beam divergence on the coherent transition radiation interference 
pattern is depicted in Figures 37-43. For all the figures, the angular information 
depicted in the X—Y plane extends to five times the predicted angle of maximum 
intensity in radians for the single foil case at the same energy. Intensity is 
measured along the Z—axis. The top of the Z—axis corresponds to the maximum 
intensity in units of charge squared divided by the speed of light per unit frequency 
and solid angle. 

The apparent inconsistency that appears in these plots is due to a forced 
asymmetry built into the program. When a beam of particles encounters a 
boundary at an oblique angle, the resulting transition radiation is assymetric in that 
one side of the pattern will be more intense than the other. This asymetry was 
simulated in the two—dimensional program by changing the sign of the observation 
angle on opposite sides of the axis of specular reflection. That technique was carried 
on into the three—-dimensional simulation. Obviously, a refinement of the method is 
needed, but the solution 1s not a simple one. The inconsistency could be removed 
by neglecting the asymetry introduced by the beam encountering the boundary at 
angle. 

Figures 37-40 depict the effect of beam divergence on the fringe pattern at an 
energy of 10 MeV. As the beam divergence is increased, the fringes become less 
visible. Note that the central cone becomes more shallow as well as the 
perpendicular component of intensity increases. All the angular information is lost 
when the beam divergence angle becomes sufficiently large for the perpendicular 


component of intensity to dominate the intensity distribution. However, the beam 
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Figure 37. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 1.41E—2 radians. The apparent inconsistency is due to a crude 
approximation of the asymmetry arising from the charged particles encountering the 
boundary at an oblique angle. The angular information depicted in the X—Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the same energy. Intensity is measured along the Z—axis. The top 
of the Z—axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 





Figure 38. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 2.83E—2 radians. Note the relative amplitudes of the outer fringes is much 
less than those of Figure 37. The angular information depicted in the X—Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the same energy. Intensity is measured along the Z—axis. The top 
of the Z-axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 39. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 5.66E—2 radians. Note that the outer fringes are no longer visible. 
Compare with Figures 37 and 38. Note also the general speading and smoothing of 
the intensity distribution. The angular information depicted in the X—Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the same energy. Intensity is measured along the Z—axis. The top 
of the Z—axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 40. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 8.49E—2 radians. The beam divergence has washed out all angular 
information. However a beam divergence angle of little more than 2.50E—2 radians 
accomplished the same effect for the single foil case. The angular information 
depicted in the X—Y plane extends to five times the predicted angle of maximum 
intensity in radians for the single foil case at the same energy. Intensity is 
measured along the Z—axis. The top of the Z—axis corresponds to the maximum 
intensity in units of charge squared divided by the speed of light per unit frequency 
and solid angle. 
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divergence angle at which this occurs is significantly higher than that for the single 
foil case. Figure 40 depicts the intensity distribution when the beam divergence 
angle is sufficiently large enough for this phenomenon to occur. 

Figures 41 and 42 depict the effect of beam divergence on the interference 
pattern at 40 MeV. Here again, the fringe parrern becomes less visible as the beam 
divergence is increased. Finally, Figure 43 depicts the case where the beam 
divergence is nearlv large enough for the angular information to be completely 
washed out. The beam divergence angle of 7.07E—3 radians required for this is well 
over twice that required of the 3.00E—3 radians for the single foil case at the same 
energy of 100 MeV. 

In general. beam divergence has the same effect on the components of intensity 
an for the coherent transition radiation case as for single foil transition radiation. 
As in the single foil case, the effects of beam divergence are amplified at higher 
energies. An increase in beam divergence decreases the visibility of the fringes. 
However, the effect of resonance causes coherent transition radiation to be less 


susceptible to the loss of angular information due to beam divergence. 


D. DIELECTRIC CONSTANTS 

The intensity of transition radiation emitted when a charged particle transits 
from one medium to another is dependent upon the dielectric constant of each 
media. It was shown in Chapter III that the intensity is directly proportional to the 
square of the modulus of the difference between the dielectric constants of the two 


media. As in the single foil case, a vacuum to medium transition was assumed. The 
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dielectric constants for gold Au, silver Ag, and aluminum Al ala wavelength of 


5890 A were used to study the effect on coherent transition radiation at 100 MeV. 





igure 41. Coherent transition radiation at, 40 MeV for an rms beam divergence 
angle of 7.07i—3 radiaus. The angular information depicted in the X—Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the sae energy. Intensity is measured aloug the Z-axis. The top 
of the Z-axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 42. Coherent transition radiation at 40 MeV for an rms beam divergence 
angle of 8.49E—3 radians. Note the decrease of the outer fringe amplitudes as 
compared to those in Figure 41. The angular information depicted in the X—Y 
plane extends to five times the predicted angle of maximum intensity in radians for 
the single foil case at the same energy. Intensity is measured along the Z—axis. The 
top of the Z—axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 43. Coherent transition radiation at 100 MeV for an rms beam divergence 
angle of 7.07E—3 radians. The beam divergence angle is nearly great enough to 
wash out all angular information. However, a beam divergence angle of only 
3.00E—3 radians was able to produce the same effect for the single foil case. The 
resonance effects in coherent transition radiation reduces the susceptibility to loss of 
angular information caused by beam divergence. The angular information in the 
X—-Y plane extends to five times the predicted angle of maximum intensity in 
radians for the single foil case at the same energy. Intensity 1s measured along the 
Z—axis. The top of the Z-axis corresponds to the maximum intensity in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 
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The results of applying the dielectric constants found in Table V of Chapter 
III for the parameters discussed here are presented in Table X for comparison with 
the results of the single foil case. The total intensity ITOT is presented in units of 


charge squared divided by the speed of light per unit frequency and solid angle. 


TABLE X. DIELECTRIC CONSTANT VS. INTENSITY FOR 
COHERENT TRANSITION RADIATION AT 100 MEV AND \=5890 A 


Metal (€)-€9)° ITOT 
ban 
Aluminum Al 138 ol 
Gold Au 1250 1964 
Silver Ag 179000 2956 


The intensities for the coherent transition radiation case are much greater 
than that for the single foil case. This effect is again attributable to resonance. It is 
clear that the intensity for the coherent case is indeed proportional to the square of 
the difference between the dielectric constants between thea two media. Figure 44 
depicts the coherent transition radiation interference pattern for silver at 100 MeV 
over an angular distribution of + 5/y radians from the Z~axis. The data is 
normalized so that the top of the Z—axis corresponds to the maximum value of 


intensity. 
E. FINITE BANDWIDTH 


Wartski [Ref. 5] found that variations in wavelength, and thus frequency, 


affected the visibility of the interference pattern fringes from the interferometer. 
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The pattern fringes became less visible as the bandwidth of wavelength increased. 
Wartski showed that the relation between intensity per unit solid angle and a 


bandwidth of wavelength could be expressed 


2 eee 2 
= 2e°(L}“AAsin 0 
[= |ry(ad)|" S I sa 


x{] ~ SORA eos ope (4.8) 


where AA is the bandwith of wavelength, L is the distance between foils, and p is 
the order of interference. An expression for the intensities of the maxima and 
minima may be found by substituting p=k+1/2 or p=K respectively. 

For this study, the dielectric constants of gold at wavelengths of 4410 A and 
5890 A was used. The coherent transition radiation distribution was studied at 30 
MeV for a single wavelength of 4410 A, and for a two percent (88 A) and fifteen 
percent (66 A) bandwidth centered on this wavelength. A similar study was done at 
150 MeV around a wavelength of 5890 A. It was found that an increase in 
bandwidth resulted in very slight fluctuations of intensity. There also appeared to 
be a slight shift in the fringe patterns for the 30 MeV case. 

Figures 45-48 depict the effect of two percent and fifteen percent wavelength 
bandwidths at 30 and 150 MeV. The data is normalized so that the top of the 
Z—axis corresponds to the maximum intensity. The X—Y plane is used to plot the 
angular information in terms of the predicted angle of maximum intensity, 1/7, for 
the single foil case. The data extends to +5/y on the X and Y axes. A close 


comparison of the fringe patterns at 30 MeV for the four percent bandwidth in 
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Figure 44. Coherent TR interference pattern for silver at 100 MeV and a 


wavelength of 5890 R. The patterns for silver and gold at the same energy and 
wavelength are nearly identical despite the difference in delectric permittivity. 

The reason for this is that the data in each case is plotted on a relative scale. The 
intensity is normalized so that the top of the Z-axis corresponnds to the maximim 
intensity. The data in the X—Y plane extends to five times the predicted angle of 
maximum intensity for the single foil case at 100 MeV. 
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Figure 45. Coherent TR interference pattern at 30 MeV for a four percent 


wavelength bandwidth centered at 4410 &. The data is normalized so that the top 
of the Z—axis corresponds to the maximum intensity. The X—Y plane is used to 
plot the angular information in terms of the predicted angle of maximum intensity, 
1/7, for the single foil case. The data extends to 5/7 on the X and Y axes. 
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Figure 46. Coherent TR interference pattern at 30 MeV for a fifteen percent 


wavelength bandwidth centered at 4410 A. A close comparison of the fringe 
patterns here with those in Figure 45 reveals that the fringe patterns become 
Slightly less visible as the bandwidth is increased. The data is normalized so that 
the top of the Z—axis corresponds to the maximum intensity. The X—Y plane is 
used to plot the angular information in terms of the predicted angle of maximum 
intensity, 1/7, for the single foil case. The data extends to 5/7 on the X and Y 
axes. 
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Figure 47. Coherent TR interference pattern at 150 MeV for a four percent 


wavelength bandwidth centered at 5890 A. The data is normalized so that the top 
of the Z-axis corresponds to the maximum intensity. The X—Y plane is used to 
plot the angular information in terms of the predicted angle of maximum intensity, 
1/7, for the single foil case. The data extends to +5/yon the X and Y axes. 
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Figure 48. Coherent TR interference pattern at 150 MeV for a fifteen percent 


wavelength bandwidth centered at 5890 A. This pattern is indistinguishible from 
that of the four percent bandwidth presented in Figure 47. The data is normalized 
so that the top of the Z—axis corresponds to the maximum intensity. The X—Y lane 
is used to plot the angular information in terms of the predicted angle of maximum 
intensity, 1/7, for the single foil case. The data extends to +5/y on the X and Y 
axes. 


on 


Figure 45 and the fifteen percent bandwidth in Figure 46 reveals that the outer 
fringes become slightly less visible as the bandwidth is increased. However, no such 
difference could be found for the same case at 150 MeV (Figures 47,48). 

In summary, a finite bandwidth of wavelength, and thus frequenccy, results in 
slight variations of intensity. An increase in bandwidth produces a small increase in 
intensity. Angular information such as fringe pattern is relatively unaffected except 


for large bandwidths. 
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RECOMMENDATIONS AND CONCLUSIONS 


The programs used in this study provide an effective tool for studying the 
characteristics of optical transition radiation to aid in interpretation of experiments. 
They provide an inexpensive method of obtaining a thorough understanding of 
optical transition radiation and its potential use for charged particle beam 
diagnostics. Such an understanding has already been promoted through the results 
of this study. 

The nature of the energy effects on angular distribution and polarization was 
thoroughly explored. Other factors affecting the intensity and distribution of 
transition radiation were also investigated It was shown that beam divergence 
affects the angular distribution of intensity, particularly at higher energies. 
However. the effects of beam divergence can be more sensitively investigated 
through the use of a stack of foils spaced at intervals comparable with the coherence 
length. Metals characterized by dielectric prermittivities much greater than one can 
be used to enhance the intensity of transition radiation. The finite bandwidths of 
the observed wavelength have little effect on the distribution of transition radiation 
unless the bandwidth is excessively large. 

Further improvements can be made on the programs used in this study. For 
the single foil and two—foil interferometer TR programs, the asymetry representing 
the effect of a charged particle encountering a boundary at an oblique angle requires 
refinement. A cosine dependence was discussed in which the asymmetry could 
beintroduced at one point and then "rotated" around the axis. The asymmetry 


factor could then be multiplied by the cosine of the angle of rotation. Some means 
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of labelling or indicating units of measure on the axes for the plots in the SON OF 
SURF3D program is also recommended. The 5ON OF SURF3D program should 
also be fully developed for commercial use. 

Finally, software should be developed to measure and present actual data in 
three dimensions for comparison with the theoretical data presented here. A 
program compatible with SON OF SURF3D under development for that purpose 
had to be abandoned due to the time constraints on this study. Such a program 
would be a necessary step in actually applying the optical transition radiation as a 


diagnostic for charged particle beams. 
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APPENDIX A 


USER'S MANUAL 


This appendix is intended to aid the operator in using the transition radiation 
programs TR3 and COHER3 and the three— dimensional graphics program SON OF 
SURF3D. While the content was written so as to be readily understandable to the 
inexperienced computer user, a basic understanding of MS—DOS procedures is 
assumed. This appendix is written in two parts. The two transition radiation 
programs, being very similar in use are treated as one unit. A description for the 
use of each was taken up to the point of invoking the use of the graphics program. 
The final section is written so that the user would be able to begin with a data file 
from either transition radiation program, plot it out in three dimensions, and vary 
all the plotting parameters in any combination so as to obtain the most satisfactory 


plot for the case at hand. 


A. THE TRANSITION RADIATION PROGRAMS 


The transition radiation programs were written in Fortran and compiled using 
the Ryan—McFarland fortran compiler. Both programs should be installed on the 
hard disk under the directory C:\FORT>. If not previously compiled, the 
programs may be compiled from within the C:\FORT> directory using the 
command RMFORT TR3 /N or RMFORT COHER3 /N. This command causes 
and object file to be formed which is used by the PLINK command, PLINK86 FI 
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TR3 or PLINKS86 FI COHER3 at the C:\FORT> prompt, to form an execute file. 
The programs are then ready for use by typing either TR3 or COHER3 for the 
program desired. 
1. Overview 
Program execution begins upon typing either TR3 or COHER3 at the 
C:\FORT> prompt. The user is then faced with the following messages in 


succession: 


ENTER NAME OF 3—D GRAPHICS FILE 

EXAMPLE: GRAPH.DAT 

[entry] 

ENTER NAME OF CORRESPONDING PARAMETER FILE 
EXAMPLE: PRMTRS.DAT 

[entry] 

WHAT IS THE DESIRED BEAM ENERGY IN MEV? 
PLEASE INCLUDE A DECIMAL PLACE. 

lentry] 


Once the above entries are made, a menu list of current parameters 
appears. This menu list for the COHER3 program differs a bit from that of the 
TR3 program as should be expected. The menu list for the TR3 and the COHER3 
programs are depicted in Figures 49 and 50, respectively. The user has the option 
of accepting the current default parameters or changing one, any combination, or all 
the parameters. In order to change any parameter, the user need only type the 
appropriate number, press return/enter, and follow the directions. In all cases, after 
the new parameter is entered, the parameter change menu reappears indicating the 
new parameters. Once all the parameters are to the user's satisfaction, the user 
accepts the current parameters by entering the appropriate number as indicated in 


order to continue with the program. 
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The last entry required of the user provides the user with a choice of 
polarization. The parallel component of intensity only, the perpendicular 
component of intensity only, or the sum of the two for the total, non—polarized 
intensity may be plotted. The decided choice is entered by typing the appropriate 


number and pressing return/enter. 


CHOOSE DESIRED POLARIZATION. 
ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION 


1 NO POLARIZATION—TOTAL INTENSITY 


2 PARALLEL. 
3 PERPENDICULAR 


Once the choice of polarization is entered, the program begins 
calculations that take approximately seven to ten minutes (assuming a math 
co—processor is installed). Progress is indicated by the program through displaying 
the current record number over which the calculations are being done. There are a 
total of 64 records corresponding to rows of data. The data is organized into 64 
rows (records) and 64 columns (fields) for use with the graphics program. All data 
for the graphics file is normalized after the intensities for all the data points are 
calculated. The output includes the graphics file of normalized data, a 
corresponding parameter file, a file containing the non—normalized intensity data in 
terms of charge squared divided by the speed of light (per unit frequency and solid 
angle), and a scratch file used for interim calculations. 

2. Required entries 
The program is executed by typing TR3 or COHER3 at the C:\FORT> 


prompt. The following message then appears upon the screen: 


ENTER NAME OF 3—D GRAPHICS FILE 
EXAMPLE: GRAPH.DAT 
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The required entry here is the name of the actual file that will be used 
by the SON OF SURF3D graphics program. It is suggested that some sort of 
mnemonic device be used for key characteristics of the file for ease in later 
identification. Table XI demonstrates the standard format used in conjunction with 


this thesis. 


TABLE XI. GRAPHICS FILE LABELLING SCHEME 


AXXXENNN.DAT A Either T for TR3 single—foil program or 
C for COHER3 Hiterfe rometer program. 
-DAT Extension required by graphics program. 
ENNN Indicated energy in MeV. 
XAX aig: key characteristics. 
Standard defaults. ie T3E10.DAT 
AN RMS Beam angle altered, N indicated 
variation number for that energy. ie 
CA2E30.DAT for second variant at 30 MeV. 
IF scattering foil present. ie TFE50.DAT 
DXX_ Dielectric constant altered. ie CDAI for 
dielectric constant of Aluminum. 
NG Data calculated over angular distance of 
N times 1/7 from axis. 1e TSGE70.DAT 
N Data normalized to maximum value for 
total intensity at 100 MeV. ie CNE90.DAT 
PERP Perpendicular component of intensity 
only selected for plots. ie TPERPES.DAT 
PAR Parallel component of intensity only 
selected for plots. ie CPARE15.DAT 
B Finite bandwidth, wavelengths altered 
COHER3 program only. ie CBE80.DAT 
TH — Foil thickness altered, COHER3 program 
only. ie CTHESO.DAT 
DL Distance between foils altered, COHER3 
program only. ie CDLE90.DAT. 


Combinations of key characteristics are also possible. For example, 
TF5GNE80.DAT indicates the presence of a scattering foil, calcualtions over an 


angular distance of +(5/7) from the axis, and data normalized to 100 MeV. 
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Upon entering the name of the graphics file, the operator is requested to 


enter the name of the associated parameters file. 


ENTER NAME OF CORRESPONDING PARAMETER FILE 
EXAMPLE; GRAPHP DAT 


The name of the parameter file should be similar enough to the name of 
the graphics file to avoid confusion at a later date. For this thesis a P for 
parameters was tagged on to the name of the graphics file. For example, 
C5GE60P.DAT or C5GP60.DAT for the parameter file corresponding to the 
C5GE60.DAT graphics file. 

One word of caution. The computer doesn't read beyond the first nine 
letters of a file name. The file name plus three—letter extension should never exceed 
twelve letters. The programs will run using filenames entered in excess of nine 
characters. However, the files in which the only difference in name occurs beyond 
the ninth character will not appear in the directories and the data will be lost. It 
ay occassionally become necessary to abbreviate the name of the graphics and 
parameter file to prevent this from happening. 

Once the name of the parameter file has been entered, the user is 
requested to enter the energy in MeV. Since the analytic solutions were obtained 
through the use of small angle approximations for angles represented by the inverse 
of the Lorentz factor, analysis at lower energies must be considered less exact. 
Energies at which the inverse of the Lorentz factor is greater than 0.1, that 
isenergies less than about 5 MeV should be treated as rough approximations. 


Angular information, which is important for purposes of beam diagnostics, is lost 
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above about 200 MeV without a polarizing filter. No such filter is built into these 
programs therefore their analvsis above 200 MeV is of little value for beam 
diagnostics. The energy entry should include a decimal point. 

Once the output files have been named and the desired energy entered, a 
menu for changing data parameters appears on the screen. Figures 49 and 50 depict 
the menus for the TR3 single—foil and the COHER3 interferometer programs 


respectively. 


The current values are: 
1 RMS beam angle is 0.001 radians/sqrt(2). 
2 IFOIL is 0. 

A value of one indicates that a scattering foil 
is present. A value of zero indicates that there 
is a target foil only. 
3 The dielectric constant of the medium is 

0.618 + 5.471. 
This program assumes a vacuum to medium transition. 
4 Angular measurement is over 0.05 radians. 
5 Normalization is to the maximum total 
intensity for the energy entered for this run. 
The option is to normalize to the maximum total 
intensity for an energy of 100 MeV. 
In all cases, parallel intensity will be 
normalized to its own maximum value for the 
energy entered for this run. 


6 Accept current parameters. 


Enter the number without the decimal for the value 
vou wish to change. Enter 6 for no changes. 


Figure 49. Parameter change menu for the single foil TR3 program 
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THE CURRENT PARAMETERS ARE: 


1 RMS beam angle in the X—Z plane is .000210 radians 
RMS beam angle in the Y—Z plane is .000210 radians 

2 IFOIL is 0. A value of 1 indicates 

that a scattering foil is present. A value of 

zero indicate absence of scattering foil. 

3 The optical constants for the back metalic foil 

are .618 + 5.47i. 

4 The optical constants for the front clear foil 

are 1.48 + 0.001. 

5 Angular measurement is over 0.05 radians 

6 Normalization is to the maximum total 

intensity for the enery entered for this run. 

The option is to normalize to the maximum total 

intensity for an energy of 100 MeV. 

7 Minimum wavelength is 4500 angstroms. 
Maximum wavelength is 4500 angstroms. 

8 Foil thickness is .005 centimeters. 

9 Distance between foils is 1.20 centemeters. 

10 Accept current parameters. 


Enter the number without the decimal for the value 
you wish to change. Enter 10 to accept parameters. 


Figure 50. Parameter change menu for interferometer COHER3 program 


The user has the option of accepting the current default parameters or 
changing one, any combination, or all the parameters. In order to change any 
parameter, the user must first type the appropriate number, then press 
return/enter. The user will then be directed to enter the new value. Each of these 
parameters will now be discussed in igieill 

Any parameter may be changed simply by entering the integer value 
that appears next to the parameter the user wishes to change. The first parameter 
option in both programs is RMS beam angle. Upon entering a 1 to change this 
parameter, the message 'PLEASE ENTER NEW VALUE FOR RMS BEAM 


ANGLE’ appears. The appropriate entry is entered in terms of radians. Not all the 
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particles within the beam travel parallel to each other. The particles tend to 
diverge at various angles within the beam. The rms beam angle is the mean of a 
Gaussian distribution of these angles of divergence. For the TR3 program it is 
assumed that the distribution of divergence in the X—Z plane (o,.) is the same as in 
the Y—Z plane (o.,) such that ae: In that case the total RMS beam angle of 
divergence is equal to ¥20. The value entered for the TR3 program should be in 
terms of radians/¥2. The COHER3 program does not make that assumption and 
allows separate values to be entered for o, and oe The values entered for the 
COHER3 program should be in radians. 

Option 2 in both programs allows for the inclusion of a separate 
scattering foil. The message 'ENTER NEW VALUE FOR IFOIL' appears when 
this option is selected. The only valid entries for this option are 0 and 1. As 
indicated. a value of 1 here indicates the presence of the scattering foil. The 
presence of the scattering foil slightly affects the distribution of the parallel 
component of TR intensity. It does not affect the perpendicular component of 
intensity. A value of 0 here indicates that a scattering foil is not present. 

The ability to change the dielectric constants of the mediums is offered 
in option 3 for the single foil TR3 program and by options 3 and 4 of the 
interferometer COHER3 program. Upon selecting any of these options, the user is 
first asked to enter the value of the real component, then asked to enter the 
imaginary component. In effect, this allows the user to change the material of the 
target foil. Option 4 of the COHER3 program allows the user to substitute another 


material for the clear front foil figured into the program. Values for dielectric 


103 


constants may be found in the CRC Handbook of Chemistry and Physics [Ref 9]. 
Both the real and imaginary parts of the dielectric constant should be entered as 
real numbers. 

Option 4 of the TR3 program and 5 of the COHER3 program allow the 
user to change the anglular range over which the data is calculated. Two different 
means are allowed for this option. The predicted angle of maximum intensty is that 
of the inverse of the Lorentz factor in radians. The user has the option of entering a 
multiple of this angle of entering a value for this angle in radians directly. Upon 


selection of this option, the following menu appears. 


Choose option 
1 Angle that is N*(1/Lorentz factor) 
2 Enter angle in radians manually. 


Choosing option 1 results in the message 'Enter N, for N*(1/Lorentz 
factor)'. The preferred entry is an integer. Choosing this option causes calculations 
of intensity to be over an integer multiple of the predicted angle of maximum 
intensity. This option aids the comparison of energy distributions over a wide range 
of energies. Choosing option 2 results in the message 'Enter desired angle in 
radians’. Obviously, the desired number is a real value of radians for the desired 
angle. Using this option to choose the same angle for a range of energies aids in 
examining the angular dependence of the intensity distribution on energy. 

The ability to normalize data to the maximum intensity value attained 
at 100 MeV is offered by options 5 and 6 of the TR3 and COHER3 programs 


respectively. This option allows for ease of study of relative intensities as a function 
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of energy. Choosing this option is irreversible. Once chosen, there is no way of 
going back to normalizing the data to its own maximum intensity for that energy. 
Upon entering the appropriate number to choose this option, the parameter change 
menu appears with the message ‘Normalization will be to 100 MeV.' in place of this 
Option. 

Option 7 of the COHER3 program allows the user to change the range 
of wavelengths over which calculations are made. Selection of this option results in 


the following messages: 


aay NEW MINIMUM WAVELENGTH. 
entry 
eny" NEW MAXIMUM WAVELENGTH. 
entry 


The user should enter the minimum and maximum wavelengths in 
angstroms. This option provides for the rare opportunity to study the effects of 
finite bandwidth on the distribution of TR intensity. The traditional development 
of transition radiation theory is built upon an evaluation of integrals that assumes a 
Dirac function for frequency bandwidth. However, a finite frequency bandwidth has 
always been encountered while obtaining real data which may account for the small 
differences between theoretical and actual data. 

Choosing option §& of the COHER3 program results in the message: 
ENTER NEW VALUE FOR FOIL THICKNESS IN CENTIMETERS. This is for 


the thickness of the front clear foil through which the beam must travel before 
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encountering the back metalic foil. This option provides the opportunity to study 
the effects of front foil thickness on the interference pattern between the forward 
radiation from the front foil and the backward radiation from the rear metallic foil. 
Option 9 of the COHER3 program allows the user to specify the 
distance between foils in centimeters. This option not only provides the user the 
ability to study coherence length, but is anecessity since the interference pattern is 
highly dependent on the distance between foils with respect to the coherence length. 
The convention has been to define coherence length as the distance a particle 


"see" the phase of its associated 


travelsin the time it takes the particle to 
electromagnetic wave change by one radian [Ref. 5]. This phase difference was 
chosen so that the particle's field and the radiation fields add coherently, hence the 
name "coherence length". 

The programs allow the user to alter any or all data parameters in any 
order or combination. Every time a user chooses an option other than ACCEPT 
CURRENT PARAMETERS a message appears which attempts to clearly specify 
the appropriate entry. Once that entry is made, the parameter change menu 
reflecting the new values reappears on the screen. This allows the user to review all 
the parameter values before continuing on with the program. Once satisfied with 
the current parameters, the user must select ACCEPT CURRENT PARAMETERS 
(option 6 for the TR3 program, option 10 for the COHER3 program) in order for 
the program to continue. 

Having accepted the current parameters, the user is requested to make 
one last selection. The transition radiation intensities are calculated by computing 


the parallel and perpendicular components separately. Total intensity is then the 
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sum of the two components. The user is provided the option of plotting either the 
total intensity or either component of intensity. This allows the user to study the 
effects of each parameter on the components of intensity separately. The user 
selects the desired coomponent (polarization) by entering the appropriate integer as 


specified on the following menu: 


ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION 
1 NO POLARIZATION—TOTAL INTENSITY 


2 PARALLEL 
3 PERPENDICULAR 


Once the selection is entered, the program will run its course, indicating 
progress by listing the current record number over which calculations are being 
made. recall that the data is organized into sixtyfour records each containg sixty 
four fields of data for use with the graphics psrogram. The intensities are 
normalized for the graphing purposes after all the other calculations have been 


completed. 


B. THE THREE—DIMENSIONAL GRAPHICS PROGRAM 


The three—dimensional graphics program used in conjunction with the 
transition radiation programs is SON OF SURF3D which was written by Don 
Gilbert of DogStar Software. SON OF SURF3D is an improved version of 
SURF3—D / SURF$87 that is not yet fully developed for public release at the time of 
this writing. The SURF3—D / SURF87 and the improved SON OF SURF3D are 
three dimensional plotting programs that are written in Turbo Pascal and use 
TurboHALO graphics routines (Ref. 10]. Turbo Pascal is a registered trademark of 


Borland International and TurboHALO is a trademark of Media Cybernetics, Inc. 
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and IMSI. The key improvement of SON OF SURF3D over SURF3—D that makes 
it useful for this application is that it is able to plot a three dimensional surface 
from a data file provided by the user. The data file used for this purpose must 


contain the following information in order 


Integer value for the number of y—rows 

Integer value for the number of x—columns 

Real value for the z—axis minimum 

Real value for the z—axis maximum 

Real value for the z—axis center (the intercept with the xy—plane) 
Real z values for each x,y point 


Unfortunately, Don Gilbert of Dogstar Software has indicated that 
development of SON OF SURF3D for public release has stopped. However a copy 


may be obtained for a fee by writing to him at the following address: 


DogStar Software 

PO Box 302 

Boomington, IN 47402 
Telephone (812) 333-5616 


A graphics device and a printer are required for the graphics program. A 
math coprocessor is preferred but not absolutely necessary. SON OF SURF3D 
supports several graphics devices and printers with specific device driver files. The 
appropriate graphics driver device file must be stored in the directory 
C:\T\HALO\ as HALODEF.DEV. The printer device driver file must be renamed 
HALODEF.PRN and placed under the directory C:\T\HALO\ . The files 
SSURF.COM and HALORTP.EXE must be installed under the directory 
C:\T\HALO. All the required files are found in the HALO.ARC file on the disk 


supplied by Dogstar Software. 
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1. USE 

Once installed, the program may be run be typing SSURF from within 
the C:\T\HALO directory. When ready, the program will display the menu 
depicted in Figure 51. The user may choose to plot any of the supplied pre—defined 
functions by pressing the appropriate letter such as C for the Cos Sqr function. 
Choosing any of the pre-defined functions in this way automatically places the 
cursor in the parameters section of the menu. Any or all the parameters may be 
then be altered prior to plotting. The user may choose to plot the program by 


pressing the function key F10 or he may choose another function by pressing Esc 


SOP Ueirons: 7—in, 1) 


C) Cos Sar z=cos( sqr(x)+saqr(y)) 

E) Error —ciee 

L) Wm Cos z=I1n(cos(x))—In(cos(y)) 

Y) Your Formula z=exp(—sar(x)/2—sqr(y)/2) 
Parameters 

Plot radius = 0.90 Z—axis Maximum=2.400 

Hatch Increment ste) X—axis Maximum=3.000 

Delta Increment == (5 Y—axis Maximum=3.000 

Hide Back Lines =O Z—axis zero =0.000 

Hide Axes =ay0 X—axis zero =0.000 

Do A—hatching = Yes Y—axis zero =0.000 | 

Do Y—hatching = No A rotation =3.000 

Points only = No B rotation =—50.0 | 

Virtual Screen = INO C rotation =—70.0 | 


| Esc=Exit F2=Print F6=File F10=Draw Arrows=Rotate | 


Figure 51. SSURF Main Menu 
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which will place the cursor back into the top part of the menu. Option Y) in the 
top part of the menu is the option of interest in that this option allows the user to 
specify a function or a complete pathname of a data file containing data as described 
above. Choosing this option does not place the cursor in the parameter section of 
the menu. Parameters may not be altered without first obtaining a plot upon 
selection of option Y). Therefore it is recommended that the user first select 
another option as described above, alter the parameters as desired, and press Esc to 
get back into the top portion of the menu before selecting option Y). In order to 
plot the three—dimensional surface of a data file, select option Y) and press 
return/enter. Then type in the complete pathname of the data file and 
pressreturn/enter again. After a few seconds, the word "file" will appear behind the 
pathname indicating that the program is ready to plot the data. The data may then 
be plotted by pressing the function key F10. Flashing indices on the z—axis indicate 
that the plot is finished. The plot may be rotated at will by pressing the arrow 
keys. The user may print the plot by pressing the F2 function key. The Esc key 
will return the display back to the main menu and place the cursor in the 
parameters section. 

The shape of the plot is largely dependent on the plotting parameters 
used to define it. Figure 52 depicts the parameters used for the plots displayed in 


this thesis. 
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Parameters 


Plot radius = 0,0 Z—axis Maximum=1.000 
Hatch Increment =e X—axils Maximum=40.00 
Delta Increment =) 05 Y—axis Maximum=40.00 
Hide Back Lines = Yes Z—axis zero =—.500 
Hide Axes = .e5 X—axis zero =0.000 
Do X—hatching = Yes Y—axis zero =0.000 
Do Y—hatching =e A rotation =3.000 
Points only = Ne B rotation =—70.00 
Virtual Screen = NG C rotation =—60.00 


Figure 52. Sample parameter values for plots. 


"Plot radius" refers to the size of the plot in the XY plane relative to the 
extent of the X and Y axes. Valid entries for "Plot radius: range from 0.01 to 1.00. 
"Hatch increment" refers to the spacing between the plotting lines. "Delta 
increment" refers to the increment used for interpolation between data points. 
Valid entries for both "Hatch increment" and "Delta increment" are 0.00 to 1.00. 

The "Hide Back Lines" and the "Hide Axes" options may be used to 
keep the plot relatively clean. The user has the option of selecting either x hatching 
or y hatching or both. Choosing neither produces the same result as "Points only", 
which overides the "X hatching" and "Y hatching" options. The "Virtual Screen" 
option supports the high—resolution capabilities of some printers. A device file 
named HALORVRI.DEV is required for this optioni. Valid entries for all these 
options are either "Yes" or "No". 

"Z—axis maximum" defines the maximum value of the vertical scale, 
regardless of the maximum value of the function or data. "X—axis maximum" and 


"Y—axis maximum" define the maximum horizontal limits of the plot, regardless of 
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the extent of the data. SON OF SURF is capable of plotting as much as a 100 by 
100 array of data centered at the intersection of the three axes. The valid entries 
for these options range from 0.01 to 100.00. 

"Z—axis zero", "X—axis zero", and "Y—axis zero" define the origin of 
respective axis relative to the origin. Selecting something other than zero for any of 
these options physically moves the entire plot. The range of valid entries were not 
listed in the scanty documentation. 

All plots may be rotated about the three axes at will, either by 
specifying the rotation in the parameters section of the menu or by using the arrow 
keys once a plot has been generated. "A rotation" appears to refer to rotation about 
the Z—axis, 'B rotation" about the X—axis, and "C rotation" about the Y—axis. 
Although the range of valid entries for these rotations are —360.00 to +360.00, all 


possible rotations can be specified within the range of —180.00 to +180.00. 


APPENDIX B 


PROGRAM LISTINGS 


A. SINGLE FOIL TR3 PROGRAM 


PROGRAM TR3 


c MODIFICATION OF PROGRAM TRADsS1, WHICH WAS WRITTEN BY 
ec. D. W. RULE NSWC, CODE MODIFIED SEPT. 1986 TO INCLUDE I PERP 
Cc modified sept. 1987 to make w(i)= itot=ipartiperp — 
c instead of the ratio of iperp to ipar 
(S MODIFIED BY WIL LONGSTAFF, NPGS, AUG. 1988 
cS FOR THREE-DIMENSIONAL ANALYSIS 
C 
C THIS PROGRAM ASSUMES A VACUUM TO MEDIUM TRANSITION 
G 
c CALC. INTENSITY DIST. AS FUNC. OF ANGLE OF TRANSITION RAD. 
c ASSUMING GAUSSIAN BEAM SCATTERING ANGLE 
C G = LORENTZ FACTOR, S = RMS BEAM ANGLE IN RADIANS /SQRT (2) 
C T = ANGLE OF OBSERVATION IN RADIANS 
INTEGER AC, BC,NNW, ROW, COL, XO, YO, COUNT, WC, LW,NC, NW, NI, NWW 
INTEGER ENTRY, CHOICE, POLAR 
REAL THETA (64), RAD2, RAD, SCALER, ZMIN, ZMAX, ZCEN, TMAX, AMAX, ANG 
REAL W(600),X(600),Y (600) ,2(600) ,TAU(600) , INT(600) , MAXY, MAXP, 
: c (4,600) ,dx (600) ,dy (600) , cy (640) , NMBR, PMAX, RMAX 
CHABPeGieR ANS*],ch*1,datf£il*20, SFILE*15, DFILS*15 
COMPLEX AI,EPS,RPAR, RPERP 
DATA P1I/3.141592654/,PS1/0.7853982/,8/0.001/, IFOIL/0/ 
DATA XN/0.618/, XK/5.47/, IPRINT/0/, INORM/1/, TWOFOIL/1.0/ 
DATA X0/32/, YO/32/, SCALER/620.258423/, ZMIN/0./, ZMAX/1.0/ 
DATA ZCEN/0.0/,ROW/64/,COL/64/,ANG/0.05/, RMAX/0.0/, XDUM/0./ 
DATA YMAX/0.0/, WMAX/0.0/, PMAX/0.0/,MAXY/0.0/,MAXP/0.0/ 
aoe (0 O, 1. 0) 
WRITE(6,*) 'ENTER NAME OF 3-D GRAPHICS FILE.' 
WRITE (6,*) 'EXAMPLE: GRAPH.DAT' 
WRITE (6,*)' ' 
READ (6,17) SFILE 
WRITE(6,*) 'ENTER NAME OF CORRESPONDING PARAMETER FILE.' 
WRITE (6,*) 'BXAMPLE: PRMTRS.DAT' 
WRITE(6,*)' ' 
READ (6,17) DFILE 
i? FORMAT (A15) 


OPEN (3, FILE='RADOUT.DAT', STATUS=' UNKNOWN’ ) 
OPEN (4, FILE=DFILE, STATUS=' UNKNOWN' ) 
OPEN (5,F ILE=SFILE, FORM='FORMATTED' ) 
OPEN (8, FILE='SCRATCH', STATUS='UNKNOWN' , ACCESS='DIRECT', 
+ RECL=256) 
REWIND 2 
REWIND 3 
WRITE (5,123) ROW, COL, ZMIN, ZMAX, ZCEN 
123 FORMAT (213, 3F8.4) 
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WRITE (6,*) 'WHAT IS THE DESIRED BEAM ENERGY IN MEV? ' 
WRITE (6,*) 'PLEASE INCLUDE A DECIMAL PLACE.' 
WRITE (6,*)' VALID ENTRIES ARE 5 MEV - 1000. MEV. ' 
READ (6, *) ENERGY 

IF (ENERGY.LT.5.OR.ENERGY.GT.1000.) THEN 

GOTO 1 

ENDIF 

G = 1.0 + ENERGY/0.511 

TMAX=1./G 


WRITE (6,*) 'The current values are:' 
WRITE (6,*)'1 RMS beam angle is ',s,' radians/sqrt(2).' 
WRITE (6, 7)" 2 LTEOUPeius . tore 
WRITE (6,7 A value of one indicates that a scattering foil' 
WRITE (6,%*)'is present. A value of zero indicates that there' 
WRITE (6, 7)) 1S a target fosl only.) 
WRITE (6,*)'3 The dielectric constant of the medium is' 
WRITE (6G, VMN) te ARK ace 
WRITE (6,%*) 'This program assumes a vacuum to medium transition. ' 
WRITE(6,*)'4 Angular measurement is over',ANG,' radians.' 
IF (ABS (RMAX-0.).LE.1E-6) THEN 
WRITE (6,*)'S Normalization is to the maximum total' 
WRITE (6,*) ‘intensity for the energy entered for this run.' 
WRITE(6,*)'The option is to normalize to the maximum total’ 
WRITE (6,*) ‘intensity for an energy of 100 Mev.' 
WRITE (6,*)'In all cases, parallel intensity will be' 
WRITE (6,%*) ‘normalized to its own maximum value for the' 
WRITE (6,*) 'energy entered for this run.' 
ELSE 
WRITE (6, *) 'Normalization will be to 100 Mev.' 
ENDIF 
WRITE (67290. a. 
WRITE (6,*)'6 Accept current parameters.' 
WRITE (6,0 ae 
WRITE (6, *) ‘Enter the number without the decimal for the value' 
WRITE (6,*) 'you wish to change. Enter 6 for no changes.' 
WRITE (6775 
READ (6,*) ENTRY 
IF (ABS (ENTRY-6) .LE.1E-6) THEN 
GOTO 3 
ENDIF 
IF (ABS (ENTRY-1) > LE. 1-6) 9 7HEN 
WRITE (6,*) 'PLEASE ENTER NEW VALUE FOR RMS BEAM ANGLE.' 
WRITE (G3.-) 7. 
READ (6. (ies 
GOTO 2 
BENDIE 
IF (ABS (ENTRY-2) .LE.1E-6) THEN 
WRITE (6,*) 'ENTER NEW VALUE FOR IFOIL.' 
WRITE (6,7) * 
READ (6, *) IFOIL 
GOTO 2 
ENDIF 
IF (ABS (ENTRY-3) .LE.1E-6) THEN 
WRITE (6, *) ‘Enter the real component: ' 
WRITE (672) 
READ (6,*) XN 
WRITE (6,*) ‘Enter the imaginary component without the "i".' 
WRITE(6, 2] 
READ (6, *) XK 
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Se 


GOTO 2 
ENDIE 
IF (ABS (ENTRY-4) . LE. 1E-6) THEN 
WRITE (6, *) ‘Choose option.’ 


WRITE (6, *)'L Angle that is N*(1/Lorentz. factor) ' 
WRITE (6,*)'2 Enter angle in radians manually.' 
Welter, (6,2) .' ' 
READ (6, *) OPTION 
ITP t{ABS (OPTION-1) .<LE.1E—-6) THEN 
WRITE (Co, ) Enter Ne for N*(l/borentz factor) ' 
WRITE ts, =)" * 
READ (6, *) NMBR 
ANG=NMBR/G 
SCALER=31./ (SIN (ANG) ) 
WRIBE6, ~) “SCALER IS “;SCALER 
ENDIF 
IF (ABS (OPTION-2) .LE.1E-6) THEN 
WRITE (6, *) ‘Enter desired angle in radians.' 
WRIDE tS, *)' 
READ (6, *) ANG 
SCALER=31./ (SIN (ANG) ) 
Weatiewo, ~) 'SCALER IS *, SCALER 
ENDIF 
GOTO 2 
ENDIF 
PPA e ot ENIRY—5).LE.1E-6) THEN 
RMAX=881.412109 
GOTO 2 
ENDIF 
WRITE (4,*) ‘The parameters for ',SFILE, 'were:' 
WRITE (4,%*)° ' 


WRITE (4,*)'  ENERGY= ',ENERGY,' Mev' 

WRITE(4,*)‘' LORENTZ FACTOR= ',G 

WRITE (4,*)' RMS beam angle: ',s,' radians/sqrt (2) .' 
Whiter). IFOIL: *,1ifoil 

WRITE (4,*)° ° 

VREPE4a,*)' A value of one indicates that a scattering foil' 


WRITE (4,*) 'was present. A value of zero indicates that there' 
WRITE(4,*)'was a target foil only.’ 


Were (4,*)' * 

WRITE(4,*)' The dielectric constant of the medium was’ 
WRITE (4,%*)' oie. 4 exe ge * 

Werle (4,*)" ' 

WRITE(4,*)' Angular measurement was over',ANG,' radians.' 


Welle (6,*) ‘CHOOSE DESIRED POLARIZATION. ’ 
WRITE (6, *) ‘ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION’ 
WRITE(6,*)' ' 
WRITE (6, *) '1 NO POLARIZATION--TOTAL INTENSITY ' 
WRITE (6, *)'2 PARALLEL. ' 
WRITE (6, *) '3 PERPENDICULAR. ' 
WMeITE(6,*)* ' 
READ (6, *) POLAR 
Pre eSiPOLAR-1).LE.1E-6). THEN 

WRITE (4,*) ' TOTAL POLARIZATION WAS SELECTED. ' 
ENDIF 
TeXABS (POLAR-2).LE.1E-6) THEN 

WRITE(4,*) ‘PARALLEL POLARIZATION ONLY. ' 
ENDIF 
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IF (ABS (POLAR-3) .LE.1E-6) THEN 
WRITE(4, *) 'PERPENDICULAR POLARIZATION ONLY.' 
ENDIF 
IF (ABS (IFOIL-1) .LE.1E-6) THEN 
TWOFOIL=0.5 
WRITE (6, *) 'SCATTERING FOIL ASSUMED TO BE PRESENT' 
ENDIF 
IF(S*G .GE. 1.0) THEN 
WRITE(6, *) 'WARNING: S MUST BE LESS THAN 1/GAMMA FOR' 
WRITE (6, *) "APPROXIMATIONS USED IN I _PERP.' 
ENDIF 
*x*xkTWFOIL ADJUSTS ASYMMETRIC TERM IN IPAR; WHEN THERE IS A 
**** SCATTERING FOIL, SET IFOIL=1, THUS TWFOIL=0.5 
B = SORTU. -S0e G2) 
S2= 2450. * Ste 
S3=SORT (P1I*S2) 
CS=COS (PSI) 
SN2=1-CS**2 
EPS=XN**2 - XK**2 4 2,*AI*XN*XK 
RPAR= (EPS*CS-CSOQORT (EPS-SN2) ) / (EPS*CS+CSORT (EPS-SN2) ) 
RPERP= (CS-CSORT (EPS-SN2) ) / (CS+CSORT (EPS-SNZ2) ) 
RSQ=CABS (RPAR) **2 
RRL=REAL (RPAR) 
DO.205l—i) 64 
WHLLe(6>.*)  becord—n it 
DO™=10 ed=10ea 
RAD2Z=REAL(( (1-XO) ~*2) + (ge) ee) 
RAD=SQRT (RAD2) 
THETA (J) =RAD/SCALER 
TE (J2GT 332) SHEN 
THETA (Co) = (THE TAC ie) 
ENDIF 
IF (12GT.322AND.J,.GT 22) te 
THETA (J) = (-THETA (J) ) 
ENDIF 
od) = THe bana) 
IF (ABS (I-1) .LT.1E-6) THEN 
X1=x* (1) 
ENDIF 
IF (ABS (1-64) .LT.1E-6.AND.ABS (J-64) .LT.1E-6) THEN 
X2=X (64) 
ENDIF 
T=THETA (J) 
FEF. = T**2/(12076*"2 43. ee 
LF (S.EO. 0 )eHEN 
D=(B**2/P1**2) * (RSQO*F+ (T/ABS(T) )) *RRL*SORT (F) 


ENDIF 
IF (S.GT.0.) CALL DIST(1,G, Byofe27 37 (VE Ole, RE Ae ERE eee 
IF (IPRINT..EQ. 1) “PRIDE ae 0 a ee 
ro) =D 
iy 3S .Gl. O20) Seen 
Z (J) =DD 


TAU (J) = (DD-D) / (DD+D) 
TAU (J) =ABS (TAU (J) ) 
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See wil) 1S ITOT = IPAR + IPERP 


W (J) =D+DD 
ENDIF 
Pee oes. 0 ..0)) THEN 
TAU TI) —1 20 
Z(J)=0.0 
W (J) =D 
ENDIF 
TPA) 26Gb. YMAX) TREN 
YMAX=Y (J) 
ENDIF 
IF (Z2(J) .GE.PMAX) THEN 
PMAX=Z (J) 
ENDIF 
IF (W(J) .GE.WMAX) THEN 
WMAX=W (J) 
AMAX=ABS (THETA (J) ) 
MAXY=Y (J) 
MAXP=2 (J) 
ENDIF 


Pea eS (POLAK-1).LE. 1E=-6) “THEN 
INT (J) =W (J) 

ENDIF 

If (ABS (POLAR-2) .LE.1E-6) THEN 
INT (J) =Y (J) 

ENDIF 

IF (ABS (POLAR-3) .LE.1E-6) THEN 
INT (J) =Z (J) 


ENDIF 
10 CONTINUE 
WRITE (8, REC=I1) (INT (COUNT) , COUNT=1, 64) 
20 CONTINUE 


Weate (3,300) ENERGY, S, XDUM,XDUM 
Weelte (3,300) X1,X2 
WRITE (3,300) YMAX, WMAX, XDUM 
MRITE(3,*) ROW, COL 
WRITE (4, *) "MAXIMUM TOTAL INTENSITY WAS ',WMAX 
WRITE (4, 7) MAX TOTAL INTENSITY OCCURED AT *,AMAX,* RADIANS’ 
Demet ae ) MAX INTENSITY PREDICTED POR. ~~; TMAX, ' RADIANS ' 
WRITE (4, *) ‘MAXIMUM PARALLEL INTENSITY WAS: ', YMAX 
WRITE (4,*) 'MAXIMUM PERPENDICULAR INTENSITY WAS: ',PMAX 
Weile (9. ) “PARALLEL INTENSITY AT MAXIMUM INTENSITY: *,MAXY 
VeellEt(4,-) PERP. INTENSITY AT MAXIMUM INTENSITY: ‘*,MAXP 
WRITE (6, *) "NORMALIZING THE INTENSITY VALUES TO 1.0' 
DO 250 NI=1, 64 
READ (8, REC=NI1) (INT (NW) , NW=1, 64) 
DO 200 NNW=1, 64 
Peron EO. 1) THEN 
Te(ABS (POLAR-2) -LE.1E-6) “THEN 
INT (NNW) =INT (NNW) /YMAX 
GOTO 100 
ENDIF 
PECABS (RMAXK—-0.) SLE. | ob=6) THEN 
INT (NNW) =INT (NNW) /WMAX 


ICY) 


ELSE 
INT (NNW) =INT (NNW) /RMAX 
ENDIF 
100 ENDIE 


& **x**NORMALIZES IPAR to unity WITH MAX VALUE OF IPAR, YMAX 


QO 


*kx*knormalizes itot=w(i) to unity with wmax 
o *x**normalizes iperp=dd with wmax also 
WRITE (3,300) X (NNW), Z (NNW), Y (NNW) , W (NNW) , TAU (NNW) 
a00 FORMAT (554614 275) 
200 CONTINUE 
WRITE (5,222) (INT (NWW) , NWW=1, 64) 
222 FORMAT (6468. 4) 
250 CONTINUE 
CLOSE "(2) 
CLOSE (3) 
CLOSE (4) 
CLOSE (3) 
CLOSE (8) 
END 


SUBROUTINE DIST (T,G,B,5,52, $3, IWEOIL, RPAR EE) Eelam 7 fe) 


c T = ANGLE OF OBSERVATION IN RADIANS 

Cc G = LORENTZ FACTOR, GAMMA 

& S = RMS SCATTERING ANGLE OR BEAM ANGLE 
Cc D = RESULTING DISTRIBUTION OF INTENSITY 


COMPLEX BETA, Z,W, Aly DIZ eZ RE EEE a co 

REAL IPERE, [PAR 

DATA PI /3214i522 6207 3MO7 17 

Al =) (O00 o> 

BETA. = 1..07Gi¥aa\ loa 

Z = BETA/ (1.414213562*S) 

IF (CABS (2) (2 LT = Getc4 ete 
©. IF MO = 1, W = BREC(Z)> IF Mo — 0, ee -eee 

CALL CEREMMO 4a) 
W=W*CEXP (2**2.) 


BLSE 
CALL LCERF (2, Week? 
cS ****LCERE RETURNS W=CEXP (2.2) ER BG7) ee Ole eo) (2) pec 
ENDIF 


D1Z = (G =e RETAy SS. 7 ea 
E12 = (G 49BETA/s7 2) eee 
D1l= REAL(D1Z) 


SORT (33) 


Bl = REAL (E1Z) 
D2 = (B*S*G)**2..% (419/83) 3 ee 2 ee ee 
B2 = B**2. * ( 41.°/8S3) "43. * S27 ier ee oe 


IPERP =~ (CABSARPEERE) = 20 aaer 2 


) 


IPAR = (CABS(RPAR)**2.) * B2 +ITWROUM* B72 EAL Re AR) AIMAG (aw) 


* / ( PI*SORT(P1I*S2) } 
RETURN 


END 
SUBROUTINE LCERE (U0, W) ter 
ANE) 1S COME a NE Wal ay Aa 
a*x*k*W IS THE VALUE OF CEXP(-—U**2) BREC(U)MOUIREA 
**** REQUIRES ABS( ARG(U) ) LESS THAN ARG=3*PI/4 


ae. @) 
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****TF CABS(U) LESS THAN SOQRT(38) THIS ROUTINE IS NOT VALID 


Gig 


10 


QaAAAAANND 


k***TER RETURNS AS 999 IF CONDITION ON ARG(U) [5 VIOLATED 
**** TER RETURNED AS -999 IF CABS(U) I5 TOO SMALL 
k*** TF THE SUBROUTINE CONDITIONS ARE VIOLATED RETURNS W= (IER, IER) 
COMP LEX U,W, US, USUM, UPROD 
DATA EP vigisha esd), naG/ 2 a5G2yente 1/1. 771245982020927 
LER=0 
X=REAL (VU) 
Y=AIMAG (U) 
TARG=ATANZ (Y, %) 
IF (ABS (TARG) _GT. ARS) THEN 
w= (999, 999) 
IER=999 
RETURN 
ENDIF 
IF ( CABS(U) .LT. 6.1644) THEN 
m= (2009. ,-999.-) 
IER=- 999 
RETURN 
ENDIF 
OS=2207 070 
USUM= (1.0,0.9) 
UPROD= (1.0,0-9) 
po 10 I=1, 3? 
UPROD=-UPROD* (2*1-1) /US 
USUM=USUM4 UPROD 
CONTINUE 
W=ouGUM/ (U*RTPT) 
RETURN 
END 
SUBROUTINE CERF (moO, A, W) 


SUBROUTINE CERF (MO, te sesanensesaaanesnneaen eee esses eeeeST 


FORTRAN SUBROUTINE FOR COMPLEX ERROR FUNCTION 
ZORTRAN SUBROUTINE COR comin ranaanaeanaansancecasese set eee00000™ 
MO = MODE OF OPERATION 
RA. = ARGUMENT (COMPLEX NUMBER) 

= FUNCTION (COMPLEX NUMBER) 
COMPLEX A, W 
DIMENSION AZ(2), Ozu S22) EF (2) 
DIMENSION TS(2), emcey, tiM(2Z), GEC ye DLO) CE (18) 
DATA S10N TS (2) 9 50000000E00/ ,CD (2) /2.086058S602 30700 1 
DATA CD 1 Oooo eo 0498687E-2/,CD (4) /1, 8542165332607 8 Ty 
DATA KD (5) 13 .27963479382361E-1/, CD (6) 75 .12675279912828E-1/ 
DATA ED 8) 3 ooo 80451 0SE-1/ -CD(B) /1.036950674 07 TED 
DATA SO oat 55437E00/ ,CD(10) /1,8689166221 8002570, 
DATA CO a0 eg 9523929E00/ ,CD (22) /3.22719383 73720 o oy 
DATA cD (11) 12 4631097626001 3E00/yCD (14) /5.S01 788730504 oo 
DATA CD (13) 14218 oe G663102600/ CD (16) /9. 451702080 6800 
DATA CD15) 17 oat eg 47 Bab +1/ CD (18) /2 724839372 S08! 
DATA CD17) 1 2 9 93324096E~2/ CE (2) /1 S927 BS76305 08 
DATA ce (3) /1 -48581625614499E-1/, CEU) 71. 33219670836245E-1/ 
DATA eee 1 ioed0co2e 78957 17 CB*O) 719 .18580959447535E-2/ 
DATA Ce a oe 4297 624E-2/ CE (8) /6, 402045380008 Se 
DATA CB 8 Oa ee G7B5E~2/ , CE (10) /3-395406584792955 27 
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eS Ce 


DATA CE (11)/2-05548099470193E-2/ , GE (12)7 1-07 047 40306 So0GE 
DATA CE(13) /4.5563489221421 9E-3/, CE (14) 7 1439 e eee oe Soe 
DATA CE (15) /3.07056139834171E-4/,CE (16) /3.78156541168541E-5/ 
DATA CE (17) /2.051735096161215-6/ (CE (18)7 2463560225602 74 72-3 
AZ (1) =REAL (A) 
AZ (2) =AIMAG (A) 
QO1l ZS=AZ (1) 7AZO) 4AZ (2) FA) 
52 Ol) =AZ Cl) AZ Gl) AZ (2) “AZ C2) 
SZ (2) =2..07AZ C17 AZ) 
OZ (1) =487 (1) 
O72) =a) 
SN=41.0 
LE (AZ (1) ) 002,003,003 
Uie2 5 OZ (1) =—AzZ. (1) 
OZ (2) =-AZ (2) 
SN=-1.0 
003 1F(ZS—-1.0) 014-014-004 
OO4 IF (ZS—-26..0) 00s700Ge00G 
005 IF (SAG 0064+ 3 747232) ela Oe 


ASYMPTOTIC EXPANSION 


006 SM(1)=0.0 
oM (2) =00 
OP (1) =tS2t ly (25725) 
OF (2)=25 22) 45.45) 
OM=0 .564189583547756*EXP (-SZ (1) ) 
TS (1) =OM*COS (-SZ (2) ) 
TS (2) =QM*SIN (-SZ (2) ) 
TM (1) =4 (TS (1) 4 O27 eect O72 ie 2 
TM (2) ==(TS (1) = O74 2 eto) * O72 Gb Zs 
PM=-0.5 
GO TO 008 
007 PM=PM+1.0 
TS (1) =TM(1) *OF (1) -TM(2) *OF (2) 
TS (2) =TMA1) * OF (2) TMi 2) > Ea) 
TM (1) =-PM*TS (1) 
TM (2) =-PM*TS (2) 
IF (ABS (SM (1) ) +ABS (TM (1))) -NEVABS (SMUD GeO To 8008 
IF (ABS (SM(2) ) +ABS (TM(2) ) .EO.ABS (SM(2)))GO TO 009 
008 SM(1)=SM(1)+TM(1) 
SM (2) =SM(2) +TM(2) 
IP (PM.bT.36275)GO TO.007] 
O09 °-1E (OZi1).. L100 1) 1607 Toma 
EF (1) =SM(1) 
EF (2) =SM(2) 
O10 -1F (MO2EO, 0) GO 10 20 
W=CMP LX (BF (1), EF (2) ) 
IF (SN.BQ.1.0) RETURN 
W=CMPLX (2.0-EF (1) ,-EF (2) ) 
RETURN 
011 EF(1)=SN* (1.0-EF (1) ) 
EF (2) =-SN*EF (2) 


120 


CLO) 


ree () 


012 


013 


014 


015 


016 


017 


018 


WeeMP BAER (1), EP (2) ) 
RETURN 


RATIONAL FUNCTION APPROXIMATION 


5M(1)=0.0 

oligZ) =0.0 

OM=0 .564189583547756*EXP (-SZ (1) ) 
Toe) —OM* COs (-SZ (2) ) 

Tene) —OM"~SIN(-SZ (2) ) 

Cit) —13S (1) *0Z (1) -TS (2) *@zut2) 
Obz?)—1S (1) *O2(2) +TS(2)*02,(1) 
DO 013 I=1,18 

mS (1) —52 (1) +Gbxd 1) 

eae oat 2) 

somo Cl) “1S (1) +TS(2)*TS(2) 
Tey —+CE(1)*1S(1)/SS 

Tiwe2) =-CE (1) *TS(2)/SS 
out) —SM (1) +TM(1) 

SM (2) =SM(2) +TM (2) 

CONTINUE 

Eivig Cb Vl) SM (1) -OF (2) *SM(2} 
BE (7) —Or (1) *SM (2) +OF (2) *SM(1) 
GO TO 010 


TAYLOR SERIES EXPANSION 


Tiel) = lee 379167 095513*OZ7Z (1) 
TM (2) =1.128379167095513*Q2Z (2) 
SM (1) =TM(1) 

SM (2) =TM (2) 

PM=0.0 

PM=PM+1.0 

DM—Z.0#PH1.0 

Poi) =the) * SZ (1) -TM (2) *SZ (2) 
Mon) =the) * SZ (2) #TM (2) *SZ (1) 
TH4(1)=-TS(1)/PM 
™(2)=-TS(2)/PM 

TS (1) =TM(1)/DM 

T5(2) =IM¢2) /DM 

IF (ABS (SM(1))+4ABS (TS(1)).NE.ABS(SM(1)))GO TO 016 
inborn sm 2) )4+ABS(TS(2)) .EBQ.ABS(SM(2)))GO TO O17 
SM (1) =SM(1)+4+TS (1) 

SM (2) =SM(2)+TS (2) 

GO TO 015 

IF (MO.NE.0)GO TO 018 

W=CMP LX (SN*SM(1),SN*SM(2) ) 
RETURN 

Eira 1). 0O-SN*SM(1) 

EF (2) =-SN*SM(2) 
W=CMPLX (EF (1), EF (2) ) 

RETURN 


MODIFIED ASYMPTOTIC EXPANSION 


zi 


Or A a 


Cc 


O13 USN=son 
GOTO 701, 


TTATEN 


TWO-FOIL INTERFEROMETER PROGRAM 


PROGRAM COHER3 
REVISION OF TEMZ, COHERI BY D.W. RULE 
REVISED BY WIL LONGSTAFF, NPGS AUG, 1988 
PLOTS NORMALIZED TR DATA IN THREE DIMENSIONS 
INCLUDES CHANGES MADE BY D.W. RULE FOR COHERENT CONTRIBUTION 
OF A CLEAR PRONT EOI: 


C*rx** PL W. RULE, NSWC MAR: 1987 REVISTOMMOES sie 
C***CALC. INTENSITY DIST. AS FUNC. OF ANGLEVOE TRANSITIONS. AG 


Ch CF aG).G) eae 


Oc) (2 ©) 


Ze 


CALCS. UNPOLARIZED DIST. AND THE PARALEED & (Eres ce iro. 
ASSUMING GAUSSIAN BEAM SCATTERING ANGLE FOR 
WARTSKI'S INTERFEROMETER 


G = LORENTZ FACTOR, S=SX=PROJECTED RMS BEAM ANGLE IN 
RADIANS IN X-Z PLANE 

SY=PROJECTED RMS BEAM ANGLE IN Y-Z PLANE IN RADIANS 

T = ANGLE OF OBSERVATION IN RADIANS MBAS. WRT 45 DEGS. 
INCLUDES BANDWIDTH OF WAVELENGTHS WLMIN TO WLMAX 
CHARACTER DATFIL*50,CH*1, SFILE*15, DFILE*15 

COMPLEX EPS, AI, 0, EPSIN, EPCOS, ZF, RTO, 22, 24, RPARF, RPERF 
COMPLEX RPARB, RPERB, EPSF, TCOEF, SDELT, RPARFSO 

INTEGER AC, BC, NNW, ROW, COL, XO, YO, COUNT, WC, LW, NC, NW, NI, NWW 
INTEGER POLAR, CHOICE, ENTRY,N, OPTION 

REAL THETA (64) , RAD2, RAD, SCALER, ZMIN, ZMAX, ZCEN, INT (64) 
REAL X(64),TAU(64), IPERP (64), IPAR(64) ,W(64),IPP, IPR 
REAL IPPF, IPPB, IPRF, IPRB, KOEFF,C (4, 64) ,DX (64) , DY (64) 
REAL TMAX, AMAX, ANG, MAXY, MAXP, NMBR, PMAX, THMX, RMAX, OMAX 


*** ENERGY = BEAM ENERGY IN MEV 

***XNB, XKB ARE REAL AND IMAG OPTICAL CONSTANTS OF METALIC (B)ACK FO 
***XNF, AKEF ARE OPTICAL CONSTANTS OF TRANSPARE IP). Ohiaiaw al: 

***XTL IS SEPARATION OF INTERFEROMETER IN CM 

***PST IS ANGLE OF BEAM WRT NORMAL TO THE FOIL, PSI=P1/4 


DATA P1I/3.141592654/, PSI/0.7853982/, THMX/0./,MAXY/0./MAXP/0./ 
DATA IFILE/3/, IPRINT/0/, INORM/1/, IFOIL/0/,SX/.000210/ 
DATA SY/.000120/,XL/1.2/, DFOIL/.005/WLMIN/4500./, DELT/7.79E-4/ 
DATA WLMAX/4500.0/, XNF/1.48/, XKF/0./, XNB/.618/, XKB/5.47/ 
DATA KK/3/,EE1/1.07, TWEOIL/1.0/, XO/32/7) 207 327 ow, Ga 
DATA COL/64/, SCALER/620.25841/, ZMIN/0.0/, ZMAX/1.0/, ZCEN/0./ 
DATA RMAX/0.0/,QMAX/0.0/, PMAX/0.0/, INT/64*0.0/ 
WRITE(6,*)' ENTER NAME OF GRAPHICS DATA FILE. ' 
WRITE (6, *)' EXAMPLE: GRAPH.DAT' 
WRITE (6,*)' ' 
READ (6,7) SFILE 
WRITE (6,*)' ENTER NAME OF CORRESPONDING PARAMETER FILE.' 
WRITE (6, *)' EXAMPLE: PRMTRS.DAT' 
WRITE(6,*)' ' 
READ (6,7) DFILE 
FORMAT (A15) 
OPEN (2, FILE='SCRATCH', STATUS='UNKNOWN', ACCESS='DIRECT', 
RECL=256) 
OPEN (3, FILE='TOUT.DAT', STATUS='UNKNOWN') 
OPEN (4, FILE=SFILE, FORM='FORMATTED') 
OPEN (5, FILE=DFILE, STATUS='UNKNOWN' ) 
WRITE (4,123) ROW, COL, ZMIN, ZMAX, ZCEN 
FORMAT (213, 3F8. 4) 
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WRITE (6,*) 'WHAT IS THE DESIRED BEAM ENERGY IN MEV? ' 
MEtie (6, ~*) PLEASE INCLUDE A DECIMAL PLACE. ' 
WMRiTihwi6e,*)' * 

READ (6, *) ENERGY 

Ge ENERGY/O.511 

THMX=1./G 

ANG = 31./SCALER 
WRITE (6,*) 'The current parameters are:' 


PRGte (5,*)' —' 
WRITE(6,*)'1 RMS beam angle in the X-Z plane is ',SX,'radians' 
WRPPE (6, *) ' RMS beam angle in the Y-Z plane is ',SY,'radians' 


meeenwe (6°) 2° [FOIL IS ', TIFOIL, ’..A value of 1 indicates’ 
WRITE (6,*) 'that a scattering foil is present. A value of' 
WRITE (6,*) 'zero indicates absence of scattering foil.' 
WRITE (6,*)'3 The optical constants for the back metalic foil' 
Wepte(6,*)"*are *, XNB,'+*,XKB; 'ia' 
WRITE (6,*)'4 The optical constants for the front clear foil' 
Wegere (6, *) “are *,XKNF "+", XKE,'1.' 
WRITE (6,*)'S Angular measurement is over ',ANG,' radians.' 
TF ERS (RMAX-0.) .LE.1E-6) THEN 
WRITE(6,*)°6 Normalization is to the maximum total' 
WRITE(6,%*) ‘intensity for the energy entered for this run.' 
WRITE(6,*) 'The option to normalize to the maximum total' 
WRITE(6,*) ‘intensity for an energy of 100 Mev.' 
WRITE (6,*)'In all casis, parallel intensity will be' 
WRITE (6,*) 'normalized to its own maximum value for the' 
WRITE (6,%*) 'energy entered for this run.' 


ELSE 
WRITE(6,*) 'Normalization will be to 100 Mev.' 
ENDIF 
WRITE (6,%*)'7 Minimum wavelength is ',WLMIN,'" angstroms.' 
WRITE (6, *) ' Maximum wavelength is ',WLMAX,' angstroms.' 


WRITE (6,*)'8 Foil thickness is ',DFOIL,' centimeters. ' 
WRITE (6,*)'9 Distance between foils is ',XL,' centemeters.' 
WRITE(6,*)'10 ACCEPT CURRENT PARAMETERS ' 
WRITE (6,7) * ' 
READ (6, *) CHOICE 
mene ot HOLCE—10).LE.1B-6). THEN 
GOTO 3 
ENDIE 
IF (ABS (CHOICE-1) .LE.1E-6) THEN 
WRITE (6,*) 'ENTER NEW VALUE FOR X-Z PLANE RMS BEAM ANGLE' 
WRITE(6,*)’° ' 
READ CG, *) SX 
WRITE (6,*) 'ENTER NEW VALUE FOR Y-2Z PLANE RMS BEAM ANGLE' 
ViloeieE(6,*)" ' 
Rep (6,*) SY 
GOTO 2 
ENDIF 
IF (ABS (CHOICE-2) .LE.1E-6) THEN 
WRITE (6, *) "ENTER NEW VALUE FOR IFOIL.' 
WRITE(6,*)'° ' 
RE30(6,*) IFOIL 
GOTO 2 
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ENDIF 
IF (ABS (CHOLCE-—3)"- CE Seo) ET) 
WRITE (6,*)' ENTER THE REAL VALUE FOR THE BACK FOIL CONSTANT. ' 
WRITE (67a 
READ (6, *) XNB 
WRITE (6, *) 'ENTER THE IMAGINARY VALUE FOR THE BACK FOIL' 
WRITE (6, *} "GONSTANT . "BO NOT VENGLUDE the. 
WRITE (Gy. aa 
READ (6,7) XEB 
COTORZ 
ENDIF 
IF (ABS (CHOICE >4) LE sir 6) THEN 
WRITE (6, *) 'ENTER THE REAL VALUE FOR THE FRONT FOIL CONSTANT. ' 
WRITE (CG, . ae 
READ (67 7) 9x 
WRITE (6, *) 'ENTER THE IMAGINARY VALUE FOR THE FRONT FOIL' 
WRITE (6, *) "GONSTANT.." DO NOT WHGECeheiie a. 


WRITE (6; *)> aa 
READIG 7 rn Ke 
COTCeZ 

ENDIF 


IF (ABS (CHOICES Ss) GbE. Ih So) aah) 
WRITE (6, 7%) Ciigose option. 
WRITE (6,*) "1 ‘Angle that 1s N* (17 lberencz taco.) 
WRITE(6,*)'2 Enter angle in radians manually.' 
WRITE(6, 7") 
READ(6,7)}) OF TION 
TE AABS (OCP TION] LE IE-6) the 
WRITE (6, *) “Enter Ni for N17 borenteseractor) 
WRETE Ga 
READ (6,*) NMBR 
ANG=NMBR/G 
SCALER=31./ (SIN (ANG) ) 
ENDIF 
TE (ABS (OP TUCN-7Z) Ee he or iin 
WRITE(6,*) 'Enter desired angle in radians.' 
WRITEs6. 25: ae 
READ (6, *) ANG 
SCALER=31./ (SIN (ANG) ) 
ENDIF 
Golo, Zz 
ENDIF 
TF (ABS (CHOICE [6))- Ge. IEG lier) 
RMAX=881.412109 
GOTO 2 
ENDIF 
TE (ABS (CHOTCEe LE. IB =6 ye ais 
WRITE (6, *) 'ENTER NEW MINIMUM WAVELENGTH. ' 
WRITENG2 
READ (6, *) WLMIN 
WRITE (6, *) 'ENTER NEW MAXIMUM WAVELENGTH. ' 
WRITE (6,7). “* 
READ (6,*) WLMAX 
GOTO 2 
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ENDIF 
PEGl eS (CHOICE-8) .LE.1TE-6) THEN 
VeRITE(6,*})" ENTER VALUE FOR FOIL THICKNESS IN CENTIMETERS. ‘ 


WRITE (6,*)* ' 
Ree (6, *) DFOIL 
GOTO 2 

ENDIF 


eee BS (CHOTCE=-9) .LE.1E-6) THEN 
PelLIE(6; *) "ENTER DISTANCE BEIWEEN FOILS IN CENTIMETERS. ' 


WRITE (6,*)‘' ' 
READ(6,*) XL 
GOTO 2 
ENDIF 
3 WRITE(5,*) 'The parameters for ‘,SFILE,' were:' 
BRIE (5,*)’ ' 
WRITE (5, *) ‘ENERGY= ', ENERGY,’ MeV' 
WRITE (5,*) ‘LORENTZ FACTOR= ',G 


WRITE (5,*)' RMS beam angle in the X-Z plane:',SX 

WRITE (5,%*)' RMS beam angle in the Y-Z@ plane:',SY 
TreelB(>, 7) TFOlL® *';IFOIL,” A VALUE OCF 1 INDICATES’ 

WRITE (5,*)* THAT A SCATTERING FOIL WAS PRESENT. A VALUE OF' 
Peewee, be 2ERO INDICATES ABSENCE OF SCATTERING FOIL.’ 
WEltE (5, *) BHE DIELECTRIC CONSTANTS FOR THE BACK METALIC FOIL: ' 


WRITE (og) enti ey KR Bae eae! 
Pelle (>, ") THEeDIELE@TRIC CONSTANTS FOR THE FRONT CLEAR FOIL:' 
WRLTE*( Or" i, | XE, 


WRITE(5,*)' Angular measurement was over ',ANG,' radians.’ 
DEQ ore X-O.) .CE.1E-6) THEN 
WRITE(5,*)' Normalization to energy selected. ' 


ELSE 
WRITE(5,*) 'Intensity was normalized to 100 Mev.' 
ENDIF 
WRITE (5,*)' MINIMUM WAVELENGTH: ',WLMIN 
WRITE(5, *) ' MAXIMUM WAVELENGTH: ',WLMAX 
WRITE (5,*) ‘Foil thickness: ',DFOIL,' centimeters.' 
WRITE (5,%*)' Distance between foils: ',XL,' centimeters.' 


Welle co, ~) CHOOSE THE BESTRED POLARIZATION. ' 
Velib«6, 7) ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION. ' 
Where. (G6, 7)! * 
Werte (6, 7) 1 NO POLORIZATION--TOIFAL INTENSITY. ' 

WRITE (6,*)'2 PARALLEL.’ 
WRITE (6,*)'3 PERPENDICULAR. ' 

WRITE (6,*)'° ' 

READ (6, *) POLAR 

IF (ABS (POLAR-1).LE.JE-6 ) THEN 

WRITE(5,%) 'TOTAL POLARIZATION WAS SELECTED’ 

ENDIF 

beso (POLAR-2).LE.1E-6) THEN 

WRITE (5, *) 'PARALLEL POLARIZATION ONLY' 

ENDIE 

Pete es (POLAR-3).LE.1E-6) THEN 

WRITE (5, %*) 'PERPENDICULAR POLARIZATION ONLY' 

ENDIF 

c****FF] is factor for incoh. addn. of single foil intensities 


N25 


c****from clear for] ffl=l.5 for loo ne senso 


c CALCULATING THE PROGRAM CONSTANTS 
S=SX 
VDG = sSORT “(1.0 1.07 Ci a 


Al—(0 20) 
IF (1FOTLo EO. metiEN 
TWEOLL=O75 
ENDIF 
BMIN=P1*XL/ (2.0* VDC*WLMAX*1 .0E-08) 
BMAX=BMIN * WLMAX /WLMIN 
IF ( WLMAX .EQ. WLMIN ) THEN 
BMAX=BMIN 
ENDIF 
S2=2 50 SoZ 
S3=SORP(PE SZ) 
DELB = BMAX- BMIN 
IF ( S .GT. 0.0 .AND. BMAX .EQ. BMIN) THEN 
DELB =1.0 
B=BMIN 
QO=1.0/S2-A1*2.0*B 
RTO=CSORT 10) 
Z2-=ie Bio G 
Z24=Q/G 
END IF 
CALL REFLEC (PSI, XNF, XKF, RPARFS, RPERES, RPARF, RPERF) 
CALL REFLEC (PSI, XNB, XKB, RPARBS, RPERBS, RPARB, RPERB) 
***RPARFS, RPERFS ARE PARALLEL & PERP. REFLECTIVITIES OF FRONT FOIL 
***WHICH IS TRANSPARENT, RPARBS AND RPERBS ARE FOR BACK METALIC FOIL 
***RPARF RPERF AND RPARB RPERB ARE THE CORRESPONDING REFLECTION COEFFS 
KOBFF = VDC**2.0*RPERBS*sv "22, 07e1*-2.0 
IF(S ~E0. 0. 0)) THEN 
B=BMIN 
ENDIF 
DO 20) T= leon 
WRITE (6, *) 'RECORD= ',I 
DO 10 J=1, 64 
RADZ=REAL (((1=xX0) 7 *2) 41) eee 
RAD=SORT (RAD2) 
THETA (J) =RAD/ SCALER 
IE(J. LT. 32) SHEN 
THETA()) =(-THEtA GD) 
ENDIF 
X (J) =THETA (J) 
IF (ABS (I-1) .LT.1E-6.AND.ABS(J-1) .LT.1E-6) THEN 
TMIN=X (1) 
TMAX= (-TMIN) 
WRITE (IFILE,15) ENERGY,G, SX, SY, XL 
WRITE (IFILE,15) TMIN, TMAX, DELT, WLMAX 
WRITE (IFILE,15) XNEF,XKF, XNB, XKB, WLMIN 
WRITE (IFILE, *) ROW, COL 
5 FORMAT (5 (2142397 ol xn} 
ENDIF 
T=THETA (J) 
LE (SS Gr. 0.0) Ven 


Cte @ 


126 


cS ero GhE ar Olley PATTERNGteeE (PERE) & IPRF (PAR) FROM (F)RONT FOII. 
ei Dol (17 Gy, VDC, SX, SY, S25 Neo, REARES, RPERFS, 
fist ink., UP PES LP RE ) 
c Peer oINGUE ePOLDePATTERH oe IEP ee (EEbRP) & IPRB (PAR) FROM (B) ACK FOIL 
TE. =" 1.50 
Abn Whol (1G, VDC, SX, 51,52, 55, Te, ReEARBS, REERBs,REARB, 
ee te RE) 
IF (BMAX .GT. BMIN) THEN 
CALL BAMDIN (KK, 1,G, S, $2, BMIN, BMAX, DPAR, DPERP) 
ELSE 
ChE Oa 2, Gore ne soe ho, oc, 24, DPAR,DPERP) 
END IF 
Looe) = 20s (i fie RB-RPARBo-DeAR/ (DELB *PI**2.0) )42.* 1IPRF 
PeeRE ()) =o 2502 (iri PPeR hObPr"DEERE/DELB ) +2.* 1PPE 
- NGreEACTORS@r. lo > INDIR AR © IPERF©1O SIMULATE 1ST SURFACE 
jp CONTRIBUTION WHICH IS INCOHERENT 
IPAR, IPERP, OR W IS INTENSITY/ (UNIT FREO., UNIT SOLID ANGLE) 
o> IN UES OF “(GHARGE**2))/ (VEL OF LIGHT) 
ELSE 


CV CEO 


Pa-1.0/G** 2 cee) ig 
Ds (71 /E Tae 2 
TPAR( J) =4m0 * REARBS/P]**2 ae Sort) ( Br ryt > 2 
PeeRe())=4..0*KOpee= Xe olN(B*ET) 7 ET )}**2 
ENDs DE 
W(J)=IPAR (J) tI PERF (J) 
c joer o-OC 0 Des CALCULATED: FROM WARTSKI °S. EXPRESSION 
Pree 2 BO: h ) PRINT’ *,- T,* re &, 
Peet POGLAR—-1).L5.i15-6 ) THEN 
LweieG Dp =W (J) 
ENDIF 
Pees (hOURR-2Z).LE.16-6) THEN 
INT (J) =IPAR (J) 
ENDIF 
IPE e Ome 3) .LE.1E-6) THEN 
INT (J) =IPERP (J) 
ENDIF 
10 CONTINUE 
Wiel TE (2, AoC teeny (COUNT) , COUNT=1, 64) 
eee ia Pei PERE Ww, TAU, THOR, BMAX, OUAX, PMAX, AMAX, MAXY, 
+ MAXP) 
DO 550 BC=1, 64 
Vier (lr iiGe oo) A (DC), IPERP(BC)o, TeAK (BC), W (BC); TAU (BC) 
35 PORMAI( a (E14. 7,.1X)) 
950 CONTINUE 
20 CONT INUE 
DO 41 AC=1, 64 
READ (2, REC=AC) (INT (WC) , WC=1, 64) 
Cc 7 waeORMOREPLACES w WITH IPAR + IPERP, I.B. THE UNPOLARIZED INT 
he CI NORM. .50. 1) THEN 
DO 200 NW=1, 64 
TE (ABS CP OLAR=7) bb .le=6)) THEN 
INT (NW) =INT (NW) /OMAX 
GOTO 200 
ENDIF 


INT (NW) =INT (NW) /RMAX 


200 CONTINUE 
ENDIF 
WRITE (4, 29) (INT (NW) , NW=1, 64) 
29 FORMAT (64F 8. 4) 
C WRITE (6, *) (INT (NW) , NW=1, 64) 
4l CONTINUE 
WRITE (5, *) ‘MAXIMUM TOTAL INTENSITY WAS ', RMAX 
WRITE (5, *) “MAX TOTAL INTENSITY OCCURED AT | Abii RADIANS. 


WRITE (5, *) "MAX INTENSITY PREDICTED FOR | THMx ae D1 ANS .. 
WRITE (5, *) ‘MAXIMUM PARALLEL INTENSITY WAS: ‘', QMAX 


WRITE (5, *) ‘MAXIMUM PERPENDICULAR INTENSITY WAS: *, PMAX 
WRITE (5, *)* PARALLEL INTENSITY AY MAXIMUM INTENS 1] 1 ey 
WRITE (5, *) ' PERP INTENSITY AT MAXIMUM INTENSI UY. Vive. 
STOP 

END 


SUBROUTINE DIST2(T,G,S,Q,B, S2,RTQ, 22, 24, DPAR, DPERP) 

T=ANGLE OF OBSERVATION, THETA RADIANS 

G= LORENTZ FACTOR 

S=SX= RMS BEAM DIVERGENCE ANGLE PROJECTED IN THE X-Z PLANE 

Q=COMPLEX FUNCTION OF S, INTERFEROMETER SPACING L, AND WAVELENGTH 

DPAR AND DPERP ARE THE INTERFERENCE TERMS IN THE DISTRIBUTION 
COMPLEX 0; AI We Why ZG 2D release ee 7 4 ee ANG: LE; 
4+ZD2,RTQ, EZEX, ZPERP 

DATA P1/3.141592654/7 MO 17 

AI=(0.0,1.0) 

Z==T7(S2Z* RTO} 

7Z.C=-Al* (21422) 

ZD=AI * (21-22) 

Z3=-A1*T/S2 

25=G-2.0*Z4 

Z10C=Z3—24 

7Z10D=Za4e 

ZEX=AIL*2.0*B* (1.0/G* *24+T**2) -4.0* (B*T) **2/0 

EZEX=CEXP (ZEX) 

CALL WERFC (ZC, ZD, WC, WD) 


CLEP Ien® 


2D2=0.5) SORT (Ply s2i ace (G/2. + 2100) “sue 
Ms +AG/ 2. =weroG pes WG = “2. RIO7 SOREL) ) 
DPAR = REAL(ZD2) 
ZPBERPS0.25° SOR? (bay 62) Gas een 
+ ( WD* (G-2.0*RTO*ZD) + WO *(G-270"RIC* 20) eae arae, SORT (2 in) 
DPERP = REAL(ZPERP) 
RETURN 
END 


SUBROUTINE WERFC (ZC, ZD, WC, WD) 
COMPLEX ZC,2ZD,WC,WD,U,W 
MO=1 
lER=0 
DO TOeiI—1),2 

IF(I .EQ. 1) THEN 
U=ZC 
ELSE 
U=ZD 
ENDIF 
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IF ( CABS(U) .LT. 6.1644) THEN 
CALL CERF (MO, U,W) 
**** CERF RETURNS W=ERFC(U) IF MO=1, SBE DAHLGREN MATHLIB TR84-143 
W=W*CEXP (U**2.) 
ELSE 
****LCERF RETURNS W=CEXP (U**2.)ERFC(U) FOR ABS(U) .GT. SQRT (38) 
CAG LOBRE (UW, TER) 
Te(iee Ne. 9) RINT ~, ARG(Uy) OUT OF BOUNDS? 
SIN CALL. 16° LEERE ' 
ENDIF 
IF (I .BQ. 1) THEN 
WC=W 
ELSE 
WD=W 
ENDIF 
CONTINUE 
RETURN 
END 
SUBROUTINE LCERE (U,W, IER) 
****1) IS COMPLEX INPUT VALUE 
oe Ne soetHe VALUBMSOE CEXP (U**2) EREC(U) OUTPUT 
****REQUIRES ABS( ARG(U) ) LESS THAN ARG=3*PI/4 
Sate weneso to) LESS THAN SORT(38) THIS ROUTINE IS NOT VALID 
****TER RETURNS AS 999 IF CONDITION ON ARG(U) IS VIOLATED 
**** TER RETURNED AS -999 IE CABS(U) IS TOO SMALL 
****ITF THE SUBROUTINE CONDITIONS ARE VIOLATED RETURNS W=(IER, IER) 
COMPLEX U,W,US,USUM, UPROD 
BAmemwet/2.141592654/, ARG/Z. 3562/7, RTP1I/1.7724538509055/ 
IBR=0 
X=REBAL (U) 
Y=AIMAG (U) 
TARG= NZ (Y,X) 
IF (ABS(TARG) .GT. ARG) THEN 
w= (999, 999) 
IER=999 
RETURN 


i 


ENDIF 

ME teers (UU) . LT. 6.1644) THEN 
WS a OO Da) 
IER=-999 
RETURN 

ENDIF 

N52 ,0*U"U0 

Usui 0,-0.0)} 

UPROD=(1.0,0.0) 

DOmrOr I=15.3/ 
UPROD=-UPROD* (2* 1-1) /US 
USUM=USUM+tUPROD 

CONTINUE 

W=USUM/ (U* RTPI) 

RETURN 

END 

SUe, OUT INE REELEC (PSI, XN, AK, REARS, RPERS, RPAR, RPER)} 

COMPLEX AI, EPS, RPAR, RPER 


NZ9 


on ene CA? Tee 


AT=(0.0, 1.0) 
CS-€OStPo) 
SNZ=1-@S4 72. 
EPS=XN**2, - XK**2. + 2.*AI*XN*XK 
RPAR= (EPS*CS-CSQRT (EPS~SN2) ) / (EPS*CS+CSQRT (EPS-SN2) ) 
RPER= (CS-CSQRT (EPS-SN2) ) / (CS+CSQRT (EPS-SN2) ) 
RPARS=CABS (RPAR) **2. 
RPERS=CABS (RPER) **2. 
RETURN 
END 
SUBROUTINE DIST(T,G,B,SX,SY,S2,S83, TWEOIL, RPARS, RPERS, 
a RPAR, IPERP, IPAR) 
T = ANGLE OF OBSERVATION IN RADIANS 
G = LORENTZ FACTOR, GAMMA 
S=SX=RMS SCATTERING ANGLE OR BEAM ANGLE PROJECTED IN X-Z PLANE 
SY=RMS SCATTERING ANGLE OR BEAM ANGLE PROJECTED IN Y-Z PLANE 
IPERP AND IPAR ARE RESULTING DISTRIBUTION OF INTENSITIES 
B = VDC IN MAIN PROGRAM, VDC = V/C, 1.E. VEL. 
+ OVER SPEED OF LIGHT 
COMPLEX BETA, Z,W,AI,D1Z,E1Z, RPERP, RPAR 
REAL IPERP, IPAR 
DATA Pl /3-1415926547, Moy 1 
Al = (O2Gyi 0} 
S=SX 
BE lAw@= 507 G — Ale i 
Z = BETA/ (1241421356275) 
IF ( CABS(Z) .LT. 6.1644) THEN 
IF MO = 1, W = ERFC(Z), IE MO = S02 We=sGer (Zz) 
CALL CERE (MO; 21) 
W=W*CEXP (Z**2.) 
ELSE 
CALL LGERE (Z2,W, 1ER) 
****LCERE ETURNS W=CEXP (Z**2) * BREC(Z) FOR ABS(Z) .GT. SORT (38} 


ENDIF 

DIZ = (G6 4= PETA SG. * 2.12 59) 

bide (Got BEINYS*"*2.)* WwW 

Dl= REAL (D1Z) 

El = REAL (E1Z) 

D2 = (B*SY*G)**2. *°( (Cl / 63) ee ee 
+ Viet (Pl asoe ~2 5) 0) 

E2 = B**2. * ( (1.783) ** 3" S27 cee Ce eS 

IPERP = REBRS = D2 7 IVECIE 

IPAR = RPARS * B2 / TWEOIL +4 B** 2. *REAL (RPAR)MAIMAG (HW) 
= / { PI*XSORT(E hss} 


***TEF SCATTERING FOIL IS PRESENT TWEOIL =90.5 AND SiIMMETRIC TERE 
***IN IPERP AND IPAR IS DOUBLED 
RETURN 
END 
SUBROUTINE NORM(T, IPAR, IPERP,W, TAU, INORM, RMAX, QMAX, PMAX, AMAX, MAXY, 
+ MAXP) 
REAL INT (64), IPAR(64), IPERP (64),W(64),TAU(64) , PMAX 
REAL QMAX, AMAX, MAXY, MAXP, RMAX 
DO 100 NNW=1, 64 
IF (IPERP (NNW) .EQ. 0.0) THEN 
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TAU (NNW) =1.0 
ELSE 
TAU (NNW) = (IPERP (NNW) -IPAR (NNW) ) / (IPAR (NNW) + IPERP (NNW) ) 
TAU (NNW) =ABS (TAU (NNW) ) 
ENDIF 
***k* ADD IPAR & IPERP FOR UNPOLARIZED INTENSITY 
W(NNW)=IPAR(NNW) + IPERP (NNW) 
IF (IPAR(NNW) .GE.QMAX) THEN 
QMAX=IPAR (NNW) 
ENDIF 
IF (IPERP (NNW) .GE.PMAX) THEN 
PMAX=IPERP (NNW) 
ENDIF 
IF (W(NNW) .GE.RMAX) THEN 
RMAX=W (NNW) 
AMAX=T 
MAXY=IPAR (NNW) 
MAXP=IPERP (NNW) 
ENDIF 
CONTINUE 
RETURN 
END 
SUBROUTINE BANDIN(KK,T,G,S,S2, BMIN, BMAX, DPAR, DPERP) 
EXTERNAL F 
DIMENSION P (5) 
DATA NN/3/,FGI2/0.0/ 
***P(1) WILL CONTAIN VALUE OF DPERP IN F(B,P) CALLS, WHILE F ITSELF 
***71S VALUE OF DPAR 
***NN=NO. OF PTS. PER SUBINTERVAL FOR GAUSSIAN INTEGRATION 
***KK=NO. OF SUBINTERVALS IN BMIN TO BMAX, INTEG. RANGE 
***F 1S NAME OF FUNCTION CALLED BY FGI : F=F(X,P) 
P(2)=T 
P (3) =G 
P(4)=S 
P5)=S2 
DPAR=FGI (BMIN, BMAX, NN, KK, F,P,FGI2) 
DPERP=FGI2 
RETURN 
END 
FUNCTION FGI(A,B,NN,K,F,P,FGI2) 
5.01.05 FGI - FORTRAN GAUSSIAN INTEGRATION ROUTINE 
5.0).05 DECK DATES AS 1S 
5.01.05 TEST CASE 36 


DIMENSION V(165),V1(90) ,W(165) ,W1 (90) ,SUM(32),P (5) , SUM2 (32) 
DIMENSION V2 (75) ,W2(75) 
EVUTVALENCE (V(1)4V2Z(1)), (V(76),Vi(1)); 
(W(1),W2(1)), (W(76),W1(1)) 
DATA. V2/ 
Ue? 14596669, 0. 7 O. 774596609, -O. 861136312, -0.339981044, 
Poe e044, 661136312,—-.906179846, -.538469310, .0, 
0.538469310, 0.906179846, -0. 932469514, -0.661209386, -0.238619186, 
Oe coer )oGo, 0.o6l209386, 0.932469514,—-0.949107912,-0.741531186, 
=Uea0aocgalS1, O. ao. 405Georosl, 0. 7491531186, 0.949107912, 


14 


+ -0.960289856, -0.796666477, -0.525532410) 0 163434672, 


0.183434642, 


+t 0.525532410, 0.796666477, 0.960289856,-07 966160240) -Uaesc0e ro 
+ -0., 6133914337 -07324253423,--0. ,. 0.324 253423 . O63 3 yi4eee 
+ 0.836031107, 0. 9681602490, -0.9739065297-0- S650 Gee, -0nc) O10 95ce, 
+t -0.433395394,-0.148874339, 0.14881433397 07 4ee ae 4 eG O40 Sore 
+ 0.865063367, 0.973906529, -0. 9782286537 -0..0 07 052000, Ue 301 5200c, 
+ —-0.519096129, -0 .269543156,40. , 0. 269548156, 02519096129; 
+ Q0.730152001, 0.887062600, 0.978228658, -0.981560634, -0.904117256, 
+ -~0.769902674, -0. 587317954, -0. 367831999, -07 25229400, 02125220407, 
+ 0.369631 499S 0 Sa 1) 354. 769902674, 02904256. 0298156052 
DATA V1/ 
+ -0. 984183055, -0. 917598300, -0. 80157000 ly 0MGa2 34s, -0 440402) 51 
+ =-0. 2304538316 Ge -230458316, 0. 44845927 51, 0.642349335- 
+ 0.8015768090) wOe or 5934007 0. 984183055, -0- 3986260505, -_0. 9284 34c60- 
+ -0,.827201315, -0.687292905, -0-5152466367 07 3 ee co 108054949, 
+ 0.108054949, 0.3191123639, 0.515248636 70. 6672 2705, 0. e220 re 
+ 0.928434884, 0.986283809, -0. 9879925135 -0 493727 3392, -0). o4e2z0coc 3. 
+ -0.724417731, -0.5710972173, -0- 3394)o tas ea 404: 
+ 0.201194094, 0.394151347, 0-5709/217 37 04 244s 06 baz eo 
+ 0.937273392, 0.987992518, -0. 989400935, -—0. 9445 ot 06 beeen 
+ -0.755404408,-0.617876244, -0.458016/778,-0- 201602). OCS 
t 0.095012500, 0.281603551, 0.458016778, O61 cee 0275240 Bee 
+t 0.865631202, 0.944575023, 0.989400935, 0.991 2630c2,-0. [eo ol loi 
+ ~0.964762256, -0.934906076, -0.896321156, -0-34920/7614, -0. 7 9440370G, 
1 -0.732182119,-0.663044267, —-0.587715757,-0 - S0Gee 07, 0-42 ssl ce 
+ -0.331868602, -0.239287362,-0.144471962, -0.048307666, 0.048307666, 
+ 0.144471962, 0.239287362, 0.331868602, 0242 et 7c) 02 50G. 770 
+ 0.587715757, 0.663044267, 0.732162 9, 04 44Gb 6) ue 56 ole 
+ 0.896321156, 0.934906076, 0.964762256, O0-9@56RI512, 905577263562, 
DATA W2/ 
+ 0.555555556, 0.88888888 9790 555555550 ee 40 oo 464 0c ol oo 
+ 0.652145155, 0.347854645. 0. 236926805) 0. 476620070 Ue ooo oe oe 
+ 0.4786286/70, 0.236926885, 0.171324492, 0.3607 Gers, 0 - 4605-3 
+ 0.467913935, 0.360761573, 0.171324492, 0.129%8496s, 0°27 52795 
+ Q0.381830051, 0.417959184, O.38183005l, 0.2737 Gss0).) 022 ace 
+ 0.101228536, 0.222381034, 9.313706646, 053676 0 3620 oe 
4+ 0.313706646, 0.222381034, 0.101220536, O20 eee, OC) cecal oa 
+ 0.260610696, 0.312347077, 0.330239355, 02322) 0. 2Z60CL0G 6 
+ 0.180648161, 0.081274388, 0.066671344, 0.149454 0-21 0ccscs 
+ 0.269266/719, 0.295524225, 0.29552422357590 5 2609266 ee eel 6G5 oo. 
+ 0.149451349, 0.066671344, 0705506856770 e250 et oe Oeil 
+ 0.233193765, 0.262804545, 0227292508 Fe0e 26 7 ee eo 
+ Q0.186290211, 0.125580369, 0 .05566856775 0-04 3336 00> 3o 
+ 0.160078329, 0:.203167427, 0.233492537. 002240 OG ee oe 
+ 0.233492537, 0.203167427, 0.1600783297 702 1G 2 00 eso 
DATA Wi/ 
+ 0.040484005, 0.092121500, 0.138873510 770-17 ee. 20 7 ole04s, 
+ 0.226283180, 0.232551553, OF 226283180 0e ZU Leela oe 
+ 0.138873510, 0.092121500, 0.04048400570. G250 40) 000s coe 
+ 0.121518571, 0.15720316/, 0.185538397,) 0.2058 eee 2 ioe, 
+ 0.215263853, 0.205198464,. 0 5185538390 7072 eee oo iy 
+ 0.080158087, 0.035119460> 0-030 7532742 0 0 0s Gee re es 220, 
+ 0.1395/0678, 0.16626920G, 0. 186leTou0 Ore Aeaee eee eo az4dZ, 
+ 0.198431485, 0.186161000, 0.-166269206,) Ui oo ie 0 1270 
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30 
40 


10 


20 


25 


+ ere eee eH + + 


Oe0703660047, 0.030753242, 0.027152459, 
DeeeG2o29 1, “0. 149595989, 02169156519, 
ero osS0GlO, 07 182603415, 0, 169156519, 
UedeolsoolZ, 0.062253524, 0.027152459, 
OPGz5572065, 0.034273863, 0.0428358698, 
Wmeaadee223, .0.072345794, 0.078193896, 
OPO 1175879, 0 7098844399, 0-095638720, 


ORUI25636720, 07099844399, 0.091173879, 

ree 5626, 02072345794, 0.065822223, 

URU4azZo 2895, YU. 054273863, 0.025392065, 
N=MINO (MAXO (NN, 3), 17) 


M=N 

IF(N.GE.17) M=32 
NADD = 0 

IF (N.LE.3) GO TO 40 
IN = N-1l 


DO 30 I=3,1N 
NADD = NADD + I 
CONTINUE 
pO Ss I=1,M 
SUM(1)=0.0 
SUM2 (1) =0.0 
CONTINUE 
t= (B-A) /FLOAT (K) 
H2=H/2. 
AA=At+h2 
DO 20 L=1,K 
DO 10 I=1,M 
Mea DD + 1 
X = H2*vV(I1) + AA 
SURIC1T) =SUM(1) +F (xX, P) 
SUM2 (1) =SUM2 (1) +P (1) 
CONTINUE 
AA=AA+H 
CONTINUE 
SUMT=0.0 
SUMT2=0.0 
DOn 2 oe), M 
II = NADD + I 
SUM = SUMI + W(1I1) *SUM(I1) 
SUMT2=SUMT2+W (11) * SUM2 (1) 
CONTINUE 
FGI=H2* SUMT 
FGI2=H2* SUMT2 
RETURN 
END 
FUNCTION F(B,P) 
DIMENSION P(5) 
COMPLEX QO, RIQ, 22, 24,AI 
AI=(0.0,1.0) 
P(1)=0. 
T=P (2) 
G=P (3) 
S=P (4) 
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Pee 253524, 
-182603415, 
5149595989; 
-007018610, 
050998059, 
S08S35Th924, 
nv JO D400 G7, 
O) ORS IA eps 200) S18) - 
0.058684093, 
Or k62Z 74395, 


oo OC ore a4 2 


moo LOBSi2Z, 
mo 2450610; 
mlzZ462897 17, 
mOUG2 74395, 
7058684093, 
moo 7692093, 
2096540089, 
92083311924, 
O2050998059, 
070070186107 


SOO O'S OO 


OO) Giant) 


001 


002 


003 
004 


S72 pro) 
O=1..0/7S2 =Al’ 2-0-8 
RTQ=CSORT (Q) 
Z2=Al1*RTO/G 
Z4=Q0/G 
CALL DIST2(T,G,S,0,8, SZ, RIO, 22,247 DE cee DEERE 
F=DPAR 
P (1) =DPERP 
RETURN 

end 

SUBROUTINE CERF (MO, A, W) 
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FORTRAN SUBROUTINE FOR COMPLEX ERROR FUNCTION 
+ i a a a a ae a a a a> ae ab ae ae ae a ae ee ae a a ae a a ae a a a ae a a a a a a a a a a a a et 
MO = MODE OF OPERATION 
A = ARGUMENT (COMPLEX NUMBER) 

= FUNCTION (COMPLEX NUMBER) 
COMPLEX A, W 
DIMENSION AZ(2)). 04 ©). 5412)7 ce (2 
DIMENSION TS (2), SM(2), IM(2})> OF (2). Chitra 
DATA CD (1) /0.00000000000000E00/, CD(2) /2.20860585601 3476-27 
DATA CD (3) /8.298069404 95687E-2/,CD(4)/1.85421653326079E-1/ 
DATA CD (5) /3.2/7963479382361E-1/,CD(6) /5.126/5279912828E-1/7 
DATA CD(7) /7.45412958045105E-1/, CD(8)7 10369506) 4ee 7 ee 
DATA CD (9) /1.40378061255437E00/,CD(10) 7 I S689 Go2 te O07 
DATA CD(11)/2246314830523929E007 ,CD(12)7 322271 oe 
DATA CD (13) /4.21534348280013E007, CDi14)7 57 50 Ber EG 
DATA CD(15)/7.19258966683102E00/7 CD (U6)7 9451 02087 G40 eeu, 
DATA CD(17)/1.25710718314784E41/,CD(16)7 1 7246353072 oe 
DATA CBE(1)/8.157230833240966-2/,CE (2) /1.59285285253437E-1/ 
DATA CE(3)/1.48581625614499E-1/,CE(4) /1.332196/70836245E-1/ 
DATA CE(5)/1.15690392878957EB-1/7,CE(6) 79. 7858095 ssw soc or-2, 
DATA CE (7) /8.05908834297624E-2/,CE(8) /6.40204538609872E-2/ 
DATA CE (9) /4.81445242767885E-2/,CE (10) /3.335406584 73295E-2/ 
DATA CE(11)/2.05548099470193E-2/,CBE(12) /1.07847403887506E-2/ 
DATA CE (13) /4.55634892214219B-3/,CE (14) /1.43984458138925E-3/ 
DATA CE(15)/3.07056139834171E-4/, CE (16) 7327791565411 695415-5/ 
DATA CE(17)/2.05173509616121E-6/7 , GE (1 Cle 2G a 64 ee 2c oe vag ba=o7 
AZ (1) =REAL (A) 
AZ (2) =AIMAG (A) 
ZS=AC CL) TWAS Y AZIZ) Aa) 
52 (1) =AZ (1) “AZ (1) -=AZ (2) “ze 
S442) =2 70" AZ) Awe 
OZ (1) =4+AZ (1) 
OZ (2) =4AZ (2) 
SN=+1.0 
TE (AZT) } 0027-003 003 
OZ (1) =-AZ (1) 
OZ(2Z) ==) 7Z) 
SN=-1.0 
TE 2 oe) 0 4 Ode od 
TEU 4e 36 0) OUS U0 GeO 
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OO) 


oe sw 


005 


006 


007 


008 


009 


010 


O11 


012 


PE sea yy Oe wGa~52(2)*S4(2)) 014,014, 012 
ASYMPTOTIC EXPANSION 


SM(1)=0.0 

ohi(2)=0.0 

QF (1) =+SZ (1) / (ZS*ZS) 

OF (2) =-SZ (2) / (ZS*ZS) 

QOM=0 .564189583547756*EXP (-SZ (1) ) 

me 1) =OM*COS (-SZ(2)) 

TS (2) =OM*SIN (-SZ (2) ) 

mt) = (1348) *OZ (1) +TS (2) *0Z(2))/2S 
Meee to OZ (2)-TS (2) *0Z(1))/ZS 
PM=-0.5 

GO TO 008 

PM=PM+1.0 

oo — i) ~Or (1) -TM (2) *OF (2) 

Moree) =1 (1) “Or (2) +TM(2) * OF (1) 

Tire y=SPM TS (1) 

TM (2) =-PM*TS (2) 

Pe ebotatitl)) +ABS (1M(1)).NB.ABS(SM(1)))GO TO 008 
PEeUABS(SM (2) )4ABS (TM(z2)) .EO.ABS(SM(2)))GO TO 009 
Shi( 1) =SM(1)4+1TM(1) 

si 2) =SMd2Z) + TM (2) 

MCE 1. 36.5)GO TO 007 
Peete) it - 0-01) GO TO 019 

BE 1) =str(l) 

EF (2) =SM (2) 

IF (MO.E0.0)GO TO 011 
W=CMPLX (EF (1), EF (2) ) 

IF (SNaB0Q.1.0) RETURN 
WeCMEEX (2 .0-EF (1),-EF (2) ) 

RETURN 

Bei aotl (1 .0-EF (1) ) 
BEeC2y SS SN* EF (2) 
W=CMPLX (EF (1), EF (2) ) 

RETURN 


RATIONAL FUNCTION APPROXIMATION 


SM(1)=0.0 

9M(2)=0.0 

Gi 09596218 95835477567 EXP (-SZ(1) ) 
is Cl) Ome COS (-SZ (2) ) 

TS (2) =QM*SIN (-SZ (2) ) 

Pay to() OZ (1) -TS (2) *02 (2) 
Ory =1541) *OZ7 (2) +TS (2) *0Z (1) 
Der O13 I=1,18 

ie) )=S2(1)+CD (T) 

TS (2) =SZ (2) 

Sos le TS (1)+1TS (2) *TS (2) 
time —=tCh(1)*TS{1)/7SS 

Ge) =-CE (i) *TS(2)/SS 

SM(1) =SM(1)+4TM(1) 


OO 


Gy. C) 


013 


01s 


O16 


O17 


018 


Oe 


SM (2) =SM(2) +TM(2) 

CONTINUE 

BE (1) —OF (1) Sua) Ob (2 )e- Src) 
EF (2)=QF (1) *SM(2) +OF (2) *SM(1) 
GO TO L010 


TAYLOR SERIES EXPANSION 


™M (1) =17128379167095513702 (1) 
TM (2) =1.128379167095513*02z (2) 
SM(1)=TM(1) 

SM (2) =TM (2) 

PM=0 .0 

PM=PM+1.0 

DM=2 .0*PM+1.0 

TS (1)=TM(1) *SZ (1) -TM (2) *SZ (2) 
TS (2) =IM 1) S22) 4 1 esol) 
™(1)=-TS (1) /PM 
™(2)=-TS(2)/PM 

TS (1) =IM Ge eM 

TS (2)=TM(2) /DM 

IF (ABS (SM(1))+ABS(TS(1)).NE.ABS(SM(1)))GO TO 016 
IF (ABS (SM(2) )+ABS(TS(2)) .EQ.ABS(SM(2)))GO TO 017 
SM(1)=SM(1)+TS(1) 

SM (2) =SM(2) +TS (2) 

GO TOC 015 

IF (MO.NE.0)GO TO 018 

W=CMPLX (SN*SM(1),SN*SM(2) ) 
RETURN 

EF (1)=1.0-SN*SM(1) 

EF (2) =-SN*SM(2) 
W=CMPLX (EF (1),EF (2) ) 

RETURN 


MODIFIED ASYMPTOTIC EXPANSION 


SNe 
GO TO ey 
END 
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